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LARGE-TIME BEHAVIOR IN NON-SYMMETRIC 
FOKKER-PLANCK EQUATIONS 

FRANZ ACHLEITNER, ANTON ARNOLD, AND DOMINIK STURZER 


Abstract. We consider three classes of linear non-symmetric Fokker-Planck 
equations having a unique steady state and establish exponential convergence 
of solutions towards the steady state with explicit (estimates of) decay rates. 
First, “hypocoercive” Fokker-Planck equations are degenerate parabolic equa¬ 
tions such that the entropy method to study large-time behavior of solutions 
has to be modified. We review a recent modified entropy method (for non- 
symmetric Fokker-Planck equations with drift terms that are linear in the po¬ 
sition variable). Second, kinetic Fokker-Planck equations with non-quadratic 
potentials are another example of non-symmetric Fokker-Planck equations. 
Their drift term is nonlinear in the position variable. In case of potentials 
with bounded second-order derivatives, the modified entropy method allows 
to prove exponential convergence of solutions to the steady state. In this ap¬ 
plication of the modified entropy method symmetric positive definite matrices 
solving a matrix inequality are needed. We determine all such matrices achiev¬ 
ing the optimal decay rate in the modified entropy method. In this way we 
prove the optimality of previous results. Third, we discuss the spectral proper¬ 
ties of Fokker-Planck operators perturbed with convolution operators. For the 
corresponding Fokker-Planck equation we show existence and uniqueness of 
a stationary solution. Then, exponential convergence of all solutions towards 
the stationary solution is proven with an uniform rate. 


1. Introduction 

Fokker-Planck equations (FPEs) describe the deterministic evolution of the prob¬ 
ability density associated to many stochastic processes [3S] . Hence, they constitute 
an important class of models in applied mathematics and an interesting object of 
study in the analysis of PDFs. This paper is concerned with the large time analysis 
of FPEs. In particular we shall analyze non-symmetric equations (corresponding 
to irreversible stochastic processes). We shall analyze the existence of unique (nor¬ 
malized) steady states and, in particular, the convergence of the time dependent 
solutions towards it. Here, the main emphasis will be put on the derivation of ex¬ 
plicit exponential decay rates (or, at least, estimates of it). Apart from an intrinsic 
mathematical interest in such decay rates, they are even relevant for the modeling of 
industrial processes, like the fiber lay-down processes in technical textile production 
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For linear, symmetric FPEs the sharp exponential decay rate equals the spectral 
gap of the generator of the evolution. But, apart from simple examples, exact 
values or good estimates of this spectral gap are rarely available. Based on the 
work of Bakry and Emery on diffusion processes [iiin], the entropy method for 
PDEs has become an important tool to study the large time behavior of wide 
classes of parabolic equations [a Ells]. The success of this approach is mainly due 
to its robustness to nonlinear perturbations and extensions [Min]- More recently 
it was even generalized to degenerate parabolic equations miiiTiE]. 

In this paper we shall focus on the large-time behavior of three classes of linear, 
non-symmetric EPEs: In ^we shall consider non-symmetric EPEs with drift terms 
that are linear in the position variable. The recent interest in these equations 
originated actually in developing entropy methods for the subclass of degenerate 
diffusions equations, or more precisely “hypocoercive” equations. But it turned 
out that this method can be viewed more naturally for non-symmetric FPEs. The 
material of this chapter will be based on the recently developed entropy method 
from [^. We shall present a review from an updated point of view and include 
several typical examples to illustrate this new method. 

In ^we shall analyze kinetic EPEs with non-quadratic potentials. Again, they 
are non-symmetric EPEs, but with a drift term that is nonlinear in the position 
variable. This will illustrate that the entropy method from 32l can be applied also 
beyond equations with linear drift, at least in perturbative settings. The material 
of this chapter is an improvement of §7 in [5]. 

21 will be concerned with FPEs with non-local perturbations. These perturba¬ 
tions will again render the evolution generator non-symmetric in an appropriately 
weighted L^-space. But, surprisingly, a wide class of non-local perturbations does 
not modify the spectrum of the underlying (standard symmetric) EPE. Hence, we 
shall use spectral methods for the large-time analysis of such models. The material 
of this chapter is an extension of m to FPEs with diffusion and drift matrices that 
are not the identity. 

2. Hypocoercive and non-symmetric Fokker-Planck 

eqnations 

In this chapter we shall study the evolution of a function f{t,x)-, t>0,xG 
under the linear FPE of the form 

(2.1) dtf = Lf := div(DV/ + Cxf), 

and subject to the initial condition f{t = 0) = fo. Without restriction of generality 
we assume that 

/o > 0 , [ fo dx=l. 

We stipulate that solutions satisfy f{t, ■) G Hence, the divergence form of 

(j2.1() implies fit, x) dx = 1 for all f > 0. In (|2.ip . the diffusion matrix D G 
is symmetric and positive semi-definite, and C G is the drift matrix. Both 

are constant in space and time. 

An important model of this class is the kinetic FPE from plasma physics [S3I1D]. 
The time evolution of the phase space probability density /(t, x, v) is governed by: 

(2.2) dtf+ v-Vxf-^xV-Vyf = vdiYyiv f) + al^^yf ] x,'!;GR";t>0. 
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Here, the position-velocity vector (x,v) plays the role of a; G d = 2n, in m- 
V, a denote (positive) friction and diffusion parameters, respectively. V = V{x) is 
a given confinement potential for the system. Next we rewrite (lOl as 


(2.3) 


ft 


diva;^„ 


0 0 

0 (TI 


^x,vf + 


-V y 

^xV + vv J 


Here, the first matrix is a singular diffusion matrix, with the identity matrix 
I G the second term is the drift. For a quadratic potential V, the kinetic 

FPE ()2.2I) takes exactly the form of (12.11) and its analysis will be covered in ^ 
The case of non-quadratic potentials is the subject of ^ 


The goal of this chapter is first to identify (under appropriate assumptions on 
D and C) the unique normalized steady state foo{x) of (12.ip . Most of all, we 
shall then study the convergence of f{t) to /oo as t —?> oo with (possibly sharp) 
exponential rates. In view of space limitations we shall mostly present only formal 
computations, which hold rigorously for regular enough solutions. But, anyhow, 
parabolic and hypocoercive FPEs regularize instantaneously to C°° (cf. Proposition 
12.61 belowl. So, regularity is actually not an issue, with the possible exemption at 
the initial time. 

2.1. Non-symmetric Fokker-Planck equations. In this section we intro¬ 
duce the notion of (non)symmetric FPEs and the relative entropy, which will be 
our main tool to analyze the large-time behavior below. For these definitions we first 
consider FPEs with x-dependent coefficients. A symmetric Fokker-Planck equation 
is defined to be of the form 

(2.4) dtf = Lif := div (D(x)[V/ + fVA{x)]), 

with a diffusion matrix D that is locally uniformly positive definite on R'^ and 
symmetric. We assume that the sufficiently regular confinement potential A satisfies 
G L^(R‘^). Then f^o ■= is the unique normalized steady state of (12.4p . The 
normalization f^o (x) da; = 1 is imposed here by changing the additive constant 
of A, which is not prescribed by (j2.4l) . The non-degeneracy of the ground state of 
Li can easily be seen from the following computation: 

(2.5) {L,f,g)H = - [ V^(^)D(a;)V(^) /oodx, 

Js.<t Joe Joo 

with {■,-)h denoting the inner product in the weighted L^-space H := 
fj^). And the right hand side of (12.51) . with f = g, vanishes iff ///oo = 
const. The quadratic form (12.51) also shows that the operator Li is symmetric in 
H (cf. §2 of [7] for more details). 

The key feature of a symmetric EPE is the gradient form of its drift vector field. 
For d > 2 we shall now consider more general drift fields, which will make the 
evolution generator non-symmetric in H. For regular diffusion matrices D(a;), the 
following equation is called a non-symmetric Fokker-Planck equation: 

(2.6) dtf = L^f := div (D(a;)[V/ + f{VA{x) + F(a;)}]) . 

Here we assume that the additional vector field F satisfies 

(2.7) div(D(a;)F(x) /oo(x)) = 0, VxGR”^, 
such that /oo = is still the unique steady state of (12.61) . 
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In typical applications, however, the FPE is given with just one drift vector 
field that is not yet split into two summands (in contrast to (j2.6l) l. In order to 
retrieve the steady state, this field then needs to be decomposed into a gradient 
part and a divergence free part (in the above sense). This task is a generalization 
of the Helmholtz-Hodge decomposition (see §2 in [4] for a typical example). Such 
a decomposition of the vector field readily yields the following decomposition of 
the operator L 2 into its symmetric part Lg in H and its anti-symmetric part Las- 
L2 = Ls + Las with 

Lsf = div (D(a:)[V/ -b fVA{x)]) , 

Lasf = div (D(a;)F(a;) /) . 

Due to (12.7(1 we have Lgfoo = Lasfoa = 0. 

Next we give a more compact form of Lg and Lag, which of course also holds for 
Li with F = 0. 


Lemma 2.1. Let D(a;) > 0, assume condition ((2.71) . and let foe = e ^ denote 
the steady state of dZl. The symmetric/anti-symmetric decomposition of L 2 then 
satisfies: 

(2.8) Lgf = div (foe D(a;)v-/-) , 

(2.9) Lasf = div (foe R(a:)V-/-) , 

J 00 ' 

where the matrices R(a:) £ are skew-symmetric and satisfy on 

(2.10) V^(R/oo) = := (DF foef - 


Proof. (12.81) is trivial, so we only discuss (1^ . The divergence-free-condition (lO) 
on the vector field G implies that there exists a matrix function B (cc) G with 

B (cc) skew-symmetric and 


( 2 . 11 ) 


G'^ix) = V^B(a;). 


Let us briefly illustrate this statement: For d = 2, 3 (12.lip simplifies to the well 
known representations of divergence free vector fields: 


B(x) = 

and, respectively, 
B(a;) = 


0 

-b{x) \ 

G = vH, 

V-L := ( 

b{x) 

0 J ’ 



0 

-hix) 

Hx) \ 

( bl{x) 

^ 3 ( 2 :) 

0 

-bi{x) 

G = curl j h2[x) 

-hix) 

bi{x) 

0 ] 

V b^lx) 


In higher dimensions, (12.111) can be verified either with differential forms (cf. §6 in 
0, m) or by Fourier transformation. 

Next we compute 

Lasf ■-= div (g = div ((V^B)^^) = (V^B)V^ = div (bV^) , 

J OQ ^ ^ J 00 ^ Joe ^ Joe' 

where we have used the skewness of B in the last two steps. Setting R := B/^^ 
yields (12.9p . □ 
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Note that Laa in (12.911 is only a first order operator - due to the skew-symmetry 
of R(x). 


As mentioned above, the main goal of this chapter is to study the convergence 
to the equilibrium for solutions to non-symmetric FPEs. To this end, our main tool 
will be the relative entropy. We define (see §2.2 of [7] for more details): 

Definition 2.2. (a) Let J be either 1R+ or K. A scalar function ip S C(J) fl 

C‘^{J) satisfying the conditions 

(2.12) = 0, V' > 0, Ip" >0, {Vf < on J 

is called entropy generator. 

(b) Let fl G Li(R‘'), /2 G with fi dx = /2 dx = 1 and ^(x) G 

J a.e. (w.r.t. the measure / 2 (dx)). Then 

(2.13) eMh)--= [ i^(^)f2dx>0 

Jm.d 

is called an admissible relative entropy of /i with respect to /2 with gener¬ 
ating function ip. 


In this definition, the term “admissible” refers to the applicability of the en¬ 
tropy method under the assumptions (I2.12p . The most important examples are the 
logarithmic entropy ei(/i|/ 2 ), generated by 

ipi {a) = (Tln(T — (T-1-1, 

and the power law entropies ep(/i|/ 2 ) with 1 < p < 2, generated by 

V’p(cr) = o-P - 1 - p(cr - 1). 

Except for quadratic entropies e ,/,2 we shall always use J = R’*'. 

The above definition clearly shows that e^(/i|/ 2 ) = 0 iff /i = / 2 . In the sub¬ 
sequent sections we shall hence try to prove that solutions f{t) to FPEs satisfy 
e)/>(/(l)l/oo) —>■ 0 as t —>■ oo. Such a convergence in relative entropy then also 
implies -convergence, due to the Csiszdr-Kullback inequality: 

ll/i - /2|lii(R<i) < e^(/i|/2) • 

This relative entropy (w.r.t. the steady state) is a Lyapunov functional for the 
evolution. As proved in §2.4 of [7] we have: 


Lemma 2.3. Let f{t) be a solution to the non-symmetric FPE (12.61) with the 
divergence-free-condition (12.71) . Then, 

d 


e4/WI/oo) = -/ (v^)^D(x)(v^) 

' Joe ' ^ /oo ' ^ J oo ' 

=■■ < 0 , 

denotes the Fisher information (of f(t) w.r.t. foo)- 


/oo dx 


dt 

(2.14) 

where I^{f{t)\f, 

We remark that the right hand side of (|2.14l) is independent of the vector field F, 
i.e. independent of Las- So, for a fixed time t, the relative entropy and its entropy 
dissipation coincide for a non-symmetric FPE and its corresponding symmetric 
FPE. 
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2.2. Hypocoercive Fokker-Planck equations. In this section we shall de¬ 
fine hypocoercivity and give some typical examples. We start with the standard 
FPE on 

(2.15) dtf = div{Vf+ xf) =: Lsf 
with the unique normalized steady state 

(2.16) fccix) = . 

As seen from (12.5L the operator L 3 is symmetric on H ■= and dissipative, 

i.e. < 0 V/ S 'D{L 3 ). Also, —L 3 is coercive in the sense that 

(-A3/,/)i/> ll/lll^, V/G{/oo}^. 

In other words, —L 3 has a spectral gap of size 1 (since 0 and 1 are the lowest 
eigenvalues of — L 3 ) and this spectral gap determines the sharp exponential decay 
of solutions towards foo- 

(2.17) 

||e^^*/o -/ooIIh < e“*||/o -/oolU, V/oGi7 with f fo dx = 1; t> 0 . 

Equilibration occurs here as a balance between diffusion and drift in p.l5l) . 

Next we consider the FPEs from m- 

aj = div(DV/+ Cx f) = Lf. 

For a singular diffusion matrix D this equation is degenerate parabolic, and the 
operator L is not coercive in where the steady state f^o will be specified 

in (12.31 below. This non-coercivity can be seen easily from (j2.8l) , when choosing 
f{x) = c - x foo{x) with a vector c G ker D. 

In spite of this lack of coercivity, such degenerate FPEs will frequently still 
exhibit an exponential convergence to equilibrium. This motivated C. Villani to 
coin the term hypocoercivity in [41] . The following definition is very general, but 
afterwards we shall only be concerned with FPEs of type (EH). 

Definition 2.4. Let be a Hilbert space. Consider a linear operator L on H 
generating a Co-semigroup {e^*)t>o- Also, consider a (smaller) Hilbert space H 
that is continuously and densely embedded in the orthogonal complement of /C := 
ker L C H (i.e. H K.^). Then, —L is called hypocoercive on H if there exist 
constants c > 1 and A > 0 such that 

(2.18) l|e'^‘/llil<ce-^‘||/||^, V/G t>0. 

In many applications to FPEs, H is a weighted L^-space, and H a weighted 
H^-space. In (I2.18P we shall typically have a leading multiplicative constant c > 1, 
while this constant is 1 in the symmetric, non-degenerate case of (12.1711 . 

Next we shall give some typical examples of such hypocoercive equations in order 
to explain the convergence mechanism. 

Example 2.1. The kinetic FPE (12.21) is non-symmetric. With a sufficiently growing 
confinement potential V(x) it is hypocoercive, and its steady state is 

foo{x,v) = ce“-t4-+^(^)] ^ 

with some normalization constant c > 0 . □ 
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Example 2.2. Next we consider the following degenerate 2D equation of form 

m-- 


(2.19) 


dtf = div 


1 0 
0 0 


V/ + 


1 —OJ 

CO 0 



=:Uf. 


For any parameter a; G R, one easily verifies that the standard Gaussian (12.161) is 
still a steady state of (12.191) , and for w 0 it is the unique normalized steady state 
/oo- For w = 0 we have drift and diffusion in the xi-direction (as in the standard 
FPE). But in the a: 2 “direction there is no equilibration. 

The term with w in (I2.19|) constitutes the rotational part of the drift matrix C 
and the anti-symmetric part of the operator L 4 . Heuristically speaking, it mixes the 
diffusive cci-direction with the non-diffusive a; 2 -direction. Hence, for fast enough 
rotations, the sharp decay rate of solutions to (I2.19|) is the average of the decay 
rates in the Xi- and a; 2 -directions. More precisely, the drift matrix C has the 
following lower bound on the real parts of its eigenvalues: 

(2.20) /X := min{Re(A) | A G cr(C)} = i for |a;| > ^ . 

As we shall show in Section [2.41 below, this lower bound determines the sharp decay 
rate ^ towards /oo- For slower rotations, however, the decay rate approaches zero 

since min{Re(A) | A G cr(C)} = | — \J\ — 

As we shall see in the decay analysis below, the decay behavior can be under¬ 
stood quite well by considering the drift characteristics corresponding to (12.191) . 
They satisfy the ODE-system Xt = —Cx. Along a characteristic, \x{t)\'^ is mono- 
tonically non-increasing, and at points with cci ^ 0 it is even strictly monotonically 
decreasing. However, when crossing the a: 2 -axis, the characteristic is tangent to 
the level curves of \x\’^ (cf. Figure [T|). As we shall see below, this implies that the 
relative entropy (e.g. e 2 (/(t)|/oo) ) may have a vanishing time derivative at certain 
points in time. 

We now indicate a possibility to obtain a strictly (and uniformly in time) decay¬ 
ing Lyapunov functional for the evolution of (12.191) . At the level of drift character¬ 
istics it is advantageous to consider (instead of |a:(t)|) the “distorted” vector norm 
{x{t), 'Px{t)) with some appropriate symmetric, positive definite matrix P. This 
P-norm will allow to realize the optimal decay of x{t) with the rate /x defined in 
(I2.20p - uniformly in time (for details, see ()2.51l) below). This idea is the essence 
of the strategy followed in m for hypocoercive equations. 

To sum up, the essence of this example is a degenerate diffusion and a rotation 
that mixes the directions. □ 


Example 2.3. Here we consider ^ 211 ) again on R.^, with the diffusion matrix 
D = diag(l, 0) and the drift matrix C = [1 0 ; 1 1]. Note that C is a (trans¬ 
posed) Jordan block. Hence, the drift characteristics (solving xt = —Cx) are here 
degenerate spirals (cf. Figure [2^). The crucial aspect of this example is that the 
asymptotic direction of these characteristics (close to x = 0 ) is not aligned with the 
diffusive Xi-direction. This again allows for equilibration as t —>■ 00 . 

One easily verifies that the steady state is given by the non-isotropic Gaussian 

fooix) = ^ 
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X ’ = - X + y 
y’ = -x 


Drift characteristic 


Level curve of P-norm 


>, 



X 


Figure 1. Drift characteristic for the 2D Fokker-Planck equation 
(12.1911 with D = diag(l, 0), C = [1 — 1 ; 1 0]: The blue spiral is 
tangent to the level curves of |a:| (black circles) when crossing the 
X 2 -axis. The red ellipse is a level curve for the “distorted” vector 
norm yj Px(t)) . For this example the optimal metric is given 
by P = [2 — 1 ; — 1 2], cf. Lemma 12.111 for the algorithm how to 
compute P. (colors only online) 

In the labeling of the two axes x means Xi and y means X 2 - 


with a normalization constant c. The contour lines of the steady state potential 
are graphed in Figure [2b. Here, the “sharp” decay rate is given by 1 — e (where 
min{Re(A) I A G a(C)} = 1). □ 

2.3. Steady states and normalized Fokker-Planck equations. In the above 
examples we saw that, to enable convergence to an equilibrium, the drift matrix C 
has to mix the diffusive and non-diffusive directions of the linear FPE mi (pro¬ 
vided D is singular). Now we give conditions on D and C such that m has a 
unique steady state: 

Definition 2.5. The operator L from (12.11) is said to fulfill condition (A) if: 

(Al) No (nontrivial) subspace of kerD is invariant under CA'. 
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Figure 2. (a) Left: Drift characteristics and flow vector field for 
the 2D FPE (gJl) with D = diag(l, 0), C = [1 0 ; 1 1]. (b) 
Right: Level curves of the quadratic potential appearing in the 
steady state, i.e. A{x) = — ln/oo(a:). 


(A2) The matrix C is positively stable (i.e. all eigenvalues have real part greater 
than zero). 

Condition (Al) is equivalent to the hypoellipticity oi dt — L (cf. §1 of [M]!. 
Moreover, it implies regularization and strict positivity of the solution to (12.11) : 

Proposition 2.6. Let condition fAl) from Definition \2.5\ hold. and let fo G 

a) Then the unique solution of 1\) satisfies / S x R'^). 

b) If fo > 0, we have f{t, x) > 0 Vt > 0, x G 

Proof. For part (a) see page 148 of [24]. Part (b) follows from the strict positivity 
of the Green’s function pertaining to (|2.1I) (see Lemma 2.5 and Theorem 2.7 in 

0 ). □ 

Condition (A2) means that there is a confinement potential that prevents the 
solution to run off to oo. Without it, there would be no steady state. Indeed, The¬ 
orems |2)8| and |2T6| will show that condition (A) is both sufficient and necessary for 
the existence of a unique normalized steady state and exponential convergence of 
solutions towards the steady state. So, for equations of type (EU, hypoellipticity 
and confinement are equivalent to hypocoercivity. 

In the following lemma we give three equivalent characterizations of the hypoel¬ 
lipticity of L that will be used for the regularization of the propagators e^*, t > 0 
(see Theorem 12.151 below, and §2 of |3T]). 

Lemma 2.7. The following three statements are equivalent, where we use k := 
rankD G {1,..., d}: 

(i) No non-trivial subspace o/kerD is invariant under CtA. 

(a) There exist constants t G {0 ,... ,d — k} and k > 0 such that 

T 

c^T>{c'^y > kI . 

3=0 


( 2 . 21 ) 










10 


F. ACHLEITNER, A. ARNOLD, AND D. STURZER 


(in) There exists a constant r S {0, ... ,d — k} such that 

(2.22) rankp^, CD5,..., C^D5] = d. 

Proof. For the equivalence of (i) and (ii) we refer to Lemma 2.3 of [5]. 

For (iii)=^(ii) let 

E := [D2, CD^,..., C^D5] e ^dx{T+i)d ^ 

Then, 

r 

^dxd 3 E ^ C'^D(C^)'’' > 0 

j=o 

has rank d and (12.211) follows. 

For (ii)^(iii) assume we had rankE^ < d. Then, 3 0 n G with E^n = 0. 
Hence, E E^u = 0 would contradict ()2.21l) . □ 

If T is the minimal constant for which (12.211) (or, equivalently, (12.221) 1 holds, then 
L fulfills Hormander’s finite rank bracket condition of order t (see [24j . Theorem 
1.1). For the explicit decomposition of the generator from (12.11) in the Hormander 
form —L = A*A + B we refer to Proposition 5 in [41]. 

As an illustration we consider the following two hypocoercive examples of ()2.1I) . 
where d = 4, A: = 2: 

Example 2.4. Let 


/1 

0 

0 

0 ^ 


/ 

1 

0 

1 


0 

1 

0 

0 



0 

1 

0 

1 

0 

0 

0 

0 

; Cl := 


-1 

0 

0 

0 

V 0 

0 

0 

0 J 


1 

0 

-1 

0 

0 J 


Here, rank[D^, CiD^^j = 4 and hence r = 1. □ 

Example 2.5. Let 


(1 

0 

0 



/1 

0 

0 

0 \ 

0 

1 

0 

0 


0 

1 

1 

0 

0 

0 

0 

0 

; C2 := 

0 

-1 

0 

1 

V 0 

0 

0 



1 0 

0 

-1 

0/ 


Here, rank[D|, C 2 DI] = 3, but rank[D 2 , C 2 DI, C 2 DI] = 4. Hence t = 2. □ 

In many works on hypocoercive equations [mill], a more restrictive assumption 
(than condition (A)) is made, namely: “No subspace of kerD should be mapped 
into kerD by C^”, which corresponds to the requirement r = 1. Let as reconsider 
the two previous examples under this aspect. In Example 12.41 Cf maps the non- 
diffusive directions from kerDi into the diffusive directions from (kerDi)^. But in 
Example 12.51 maps the non-diffusive direction (0, 0, 0, 1)^ G kerD 2 still onto 
the non-diffusive direction (0, 0, 1, 0)^ G kerD 2 . But in a second step, we have 
C^(0, 0, 1, 0)^ = (0, 1, 0, —1)^, which has a nontrivial component in the diffusive 
subspace (kerD 2 )'‘“. 
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Next we discuss the existence of a steady state to (12. ip (for the proof cf. §1 of 
[3l] or Th. 3.1 in [5]): 

Theorem 2.8. There exists a unique steady state foo G of 112.1\} fulfilling 

fgd foo dx = 1 iff condition (A) holds. 

Moreover, this steady state is of the (non-isotropic) Gaussian form 

(2.23) foo{x) =CKex.p(^-l^—^ — 

where K is the unique, symmetric, and positive definite solution to the continuous 
Lyapunov equation (cf. pS] J 

(2.24) 2D = CK-fKC^, 

and ck = (27r)“^ (det K)“5 is the normalization constant. 

With the steady state at hand, we now give the decomposition of the operator 
L from ([24)) : 

Lemma 2.9. Let L satisfy condition (A). Then, its symmetric/anti¬ 
symmetric decomposition satisfies: 

(2.25) L,/ = div(/ooDV^) , 

(2.26) Lasf = div (/ooRV-/-) , 

Joe' 

with R := i(CK - KC^) 0. 

Proof. To reduce this result to Lemma [2dl we first compare (12.11) to (12.61) : The 
drift vector field Cx of CH) corresponds to D{VA + F}. Hence, we have with 
ffHM and VH = K-ix: 

Di^(a;) = Cx- DVH(x) = [C - i(CK -f KC^)K-i]a; 

(2.27) = i(CK-KC^)K'ia:. 

To verify the divergence-free-condition (lO) we compute 

div(DF(a:) foo{x)) 

= i Tr([CK - KC^]K-i) f^ - [CK - KC^] = 0, 

due to the skew-symmetry of CK — KC^. 

Next we verify the condition ()2.10l) : 

(V^(R/oo))^ = -RV/oo = i(CK - KC^)K-ix/oo = T>F{x)foo{x ), 

where we used (I2.27P in the last step. The claims (I2.25p . (I2.26P then follow from 
Lemma [O 

Finally we prove that R ^ 0. Otherwise (12.241) would imply D = KC^, and 
hence kerD = kerC^, which would contradict condition (A). □ 
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The result R 0 shows that hypocoercive FPEs of form (|2.1I1 are always non- 
symmetric. 

Next we shall bring the hypocoercive FPEs (ED to a normalized form, which 
will simplify our computations below. With its steady state given in (|2.23l) we 
introduce, as a first step, the coordinate transformation y := x G With 

5(2 /) ■= /(a/K?/), this transforms (12. ip to 

dtg = divy(DVy5 + Cyg), 

with D = and C = y/K ^CVK. a simple computation, using 

(12.241) shows that 

D = C,, 

where Cg := i(C + C^) denotes the symmetric part of C. Clearly, the transformed 
steady state reads goc{y) = with some normalization constant c > 0. As 

a second step we rotate the coordinate system to diagonalize the diffusion matrix: 
For an orthogonal matrix U G let D := U^DU = diag(di, ..., dk, 0,..., 0), 

where k = rankD. We set z := and h{z) := g{\Jz), which satisfies 


(2.28) 


dth = divz(DV2/i + Cz h). 


The symmetric part of the new drift matrix, C = U^CU, again satisfies D = Cg. 
Since the matrices C and C are similar, we have cr(C) = cr(C), which will be the 
quantity that determines the decay rate of a hypocoercive FPE. We also note that 
hoa{z) = ce“(^( with some normalization constant c. 

We remark that the above Examples [221 ED and l2.5l are already of this normal¬ 
ized form, but Example ES] is not. The above normalization brings Example ESI to 
the form 

' 2 0 \ ^ , (2 1 

0 0 ) ( -1 0 


dth = div2 


z h 


Scaling time by a factor i shows that this equals the FPE in Example 12.21 with the 
rotation parameter w = — i, which is a limiting case in (12.201) . 

For the rest of this chapter we shall always assume that the FPEs are normalized 
as in (12.281) . So, the matrices in (12.11) will satisfy D = Cg with D being diagonal, 
which implies K = I. 

2.4. Modified entropy method. To start with, let us very briefly review the 
standard entropy method for FPEs (cf. [101 [^ for symmetric FPEs and mm 
for non-symmetric FPEs): In a first step one establishes a differential inequality 
between the Fisher information (12.14p of a solution f{t) and its time derivative, 
which yields exponential decay of the Fisher information. We give the result for 
symmetric FPEs: 

Lemma 2.10. Let f(t) be the solution to ()2.4p with a constant dijfusion matrix 
D. Let the coefficients of this FPE satisfy the following Bakry-Emery condition for 
some Ai >0; 

d^A, 


dx"^ 


(x) > AiD , Vx G 


(2.29) 
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Also, let the initial condition satisfy /i/i(/o|/oo) < oo. Then, the Fisher information 
decays exponentially: 

(2.30) Idfit)\foo) < e-"^^*/v,(/o|/oo) , t>0. 


Proof. After a lengthy computation the time derivative of the Fisher information 
can be written as follows (for scalar diffusions D{x) cf. Lemma 2.13 in [7], and for 
the generalization to non-symmetric FPEs (12.61) cf. Lemma 2.3 in m)- Using the 
notation u := we have: 

= -2 / dx-2 [ Tr(XY)/oo dx 

at J^d \Joo-' Js^d 

(2.31) <-2Ai/ u^-Dufoodx = -2XiI^{fit)). 

JRd ^foa^ 


In the last estimate we used the Bakry-Emery condition (12.291) and Tr (XY) > 0. 
Here, the two matrices X, Y G are defined as follows: 


(2.32) 


X(x) := 


f r 

r' \ 

f f(x) \ 

V r' 

) 

V/oo(a:)/ 


> 0, Vx G 


since det X = — (i/'")^ > 0 for admissible relative entropies (cf. (j2. 121) 1. 


(2.33) 


Y(x) := 


Tr[(Dt)^] 

OX 



> 0, Vx G 


due to the Cauchy-Schwarz inequality. The differential inequality (12.311) for I,p{f{t)) 
implies (I2.30L and it can be written equivalently as e"{t) > —2Xie'{t) (with e{t) := 

e^(/(i)l/oo)). □ 


In the second step of the entropy method one proves the exponential decay of 
the relative entropy ()2.13|) of /(f) w.r.t. /oo. To this end one integrates (I2.31|) from 
f to oo, which yields the entropy inequality 

(2.34) -^e^(/(f)|/oo) <-2Aie^(/(f)|/oo), Vf>0. 

Hence, the relative entropy decays exponentially: 

(2.35) ev,(/(f)|/oo)<e-2^^‘e^(/o|/oo), Vf>0. 


Next we illustrate how the situation changes from a symmetric FPE to a non- 
symmetric or even hypocoercive FPE. In a symmetric FPE with D > 0, the rel¬ 
ative entropy is a convex function of time, and the entropy dissipation satisfies 
e^(/|/oo) < 0 for all probability densities f ^ foe (cf. Figure [3]). For a hypoco¬ 
ercive FPE with a singular diffusion matrix D, however, eft) is not convex. In 
fact, it decays in a “wavy” fashion, and it may have horizontal tangents at equally 
spaced points in time (cf. Eigure|4]). This oscillatory behavior is also known from 
space-inhomogeneous kinetic equations (cf. §3.7 of [40]; and m for a numerical 
study on the Boltzmann equation). 

So we observe that the entropy dissipation e(^(/|/oo) may vanish for certain 
probability densities / 7^ foe- This can also be seen from the form of the Fisher in¬ 
formation in ([2T4|): choose/(x) = (l-|-c-x)/oo(a:) with a vector c G kerD. Hence, 
an entropy inequality of the form (12.341) cannot hold for degenerate, hypocoercive 
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Figure 3. Prototypical behavior of the logarithmic relative en¬ 
tropy ei(/(t)|/oo) (solid red curve), its first (dotted black), and sec¬ 
ond time derivative (dashed blue) for a non-degenerate, symmetric 
FPE: The inequalities e' < —2Ae, e" > —2Ae' can be obtained, 
(colors only online) 


FPEs! We also see: While the Eisher information /i/,(/(t)|/oo) is a Lyapunov func¬ 
tional for symmetric FPEs, its non-monotonicity in the hypocoercive case makes it 
“useless” there. As a remedy, we present now a modified entropy method for FPEs 
of the form m, normalized as introduced in 112.31 

Since the above problems stem from the singularity of D, we now define a mod¬ 
ified entropy dissipation functional: 

(2.36) S^{f):= [ rii-) 

J'MA ^ Joo' ^ Joo' 

where the positive definite matrix P G still has to be determined. Note that 
the only difference to the Fisher information is the replacement of the matrix D 
there by P here. This auxiliary functional will take over the role of in the first 
step of the entropy method. So, our goal is to derive a differential inequality between 
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Figure 4. Prototypical behavior of the logarithmic relative en¬ 
tropy ei(/(t)|/oo) (solid red curve), its first (dotted black), and 
second time derivative (dashed blue) for the degenerate, hypocoer- 
cive FPE from Example 12.21 with D = diag(l, 0), C = [1 — 1; 1 0] 
: The inequalities e' < —2Ae, e" > —2Ae' are wrong, in general, 
(colors only online) 


S^{f{t)) and -^S^{f{t)) for a “good” choice of P > 0. Then, once exponential 
decay of ^^(/(t)) is obtained, the trivial estimate P > cpD (with some cp > 0) 
implies 

and also exponential decay of follows. 

The key question for using the modified entropy dissipation functional S^{f) is 
how to choose the matrix P for a given, normalized FPE. To determine P we shall 
need the following algebraic result: 

Lemma 2.11. For any fixed matrix Q S let /i := min{Re{A}|A is an eigen¬ 

value o/Q}. Let {Am|l < rn < mo} be all the eigenvalues of Q with Re{Am} = fJ,, 
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only counting their geometric multiplicity. 


(i) If all Xm, rn € {1,..., mo}, are non-defectiv^, then there exists a symmet¬ 
ric, positive definite matrix P G with 

(2.37) PQ + Q^P >2fiP. 

(ii) If Xm is defective for at least one m G {1,..., mg}, then for any e > 0 there 
exists a symmetric, positive definite matrix P = P(£) G with 

(2.38) PQ + Q^P > 2(/i - e)P . 


Proof. Here we only give the proof for the case that Q is not defective (and hence 
diagonalizable) and refer to Lemma 4.3 in for the general case. Let wi,... ,Wd 
denote the eigenvectors of Q^. Then one can choose P as a weighted sum of the 
following rank 1 matrices: 

d 

(2.39) P := bj Wj 0 , 

i=i 

with bj G j = 1,... ,d. As is a basis of P is positive definite. 

If any Wj is complex, its complex conjugate WJ is also an eigenvector of Q^, since 
Q is real. By taking the same coefficient bj for both, we obtain a real matrix P. 
Apart from this restriction, the choice of bj > 0 is arbitrary. For P from (12.3911 . we 
have 

d d 

PQ + Q^P = ^ bj {Xj + Xj )wj ® Wj'^ > 2fi ^ bj Wj (Si wj'^ = 2p,P . 

i=i i=i 


□ 


We remark that P is, in general, not unique, not even up to a multiplicative 
constant. But this will be irrelevant for the decay rate of FPEs. 


Applying this lemma to Q := C now yields exponential decay of the functional 
Sjj{f(t)), defined with the matrix P from the above lemma: 

Proposition 2.12. Assume condition (A). Let if generate an admissible entropy, 
let f be the solution to with an initial state satisfying S',^(/o) < oo, and let 

pL := min{Re{A}|A is an eigenvalue of C} (which is positive by condition (A)). Let 
{Am|l < m < mg} be the eigenvalues of C with Re{Am} = P defined 

as in Lemma \2.11\ 

(i) If all Xm, 1 < m < mg, are non-defective, then 

Sj,{fit)) < S4fo)e-^^\ t > 0. 

(ii) If Xm is defective for at least one m G mg}, then 

Sp,{f{t),e) < Sj,{fQ, t > 0, 

for any s G (0,//). Here, Stp{f,e) denotes the modified entropy dissipation 
functional (|2.36p with the matrix P = P(e). 

^An eigenvalue is defective if its geometric multiplicity is strictly less than its algebraic 
multiplicity. 
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Proof. In a tedious computation the time derivative of the functional S{'if>{f{t)) can 
be written as follows (cf. Proposition 4.5 in (5). Using the notation u := -f- we 

have: 

-^S^{f{t))=-2 [ V{t-) u^[PC + C^P]u/oo dx 

(2.40) -2 [ Tr(XYp)/oodx 

<-2k[ u^Pm/oo dx = -2KS^{f{t)), 

jR<i ^Joo' 


where k := for case (i), and k ^ — e ioi the defective case (ii). In the last 
estimate we used the matrix inequality (12.371) in case (i) and (12.381) for case (ii). This 
inequality replaces the Bakry-Emery condition (12.291) used in the standard entropy 
method (compare to the estimate (12.311) 1. In (12.401) we also used Tr(XYp) > 0, 
where the matrix X is defined in (12.321) . and the matrix Yp S is now defined 
as follows: 


Yp(a:) 


Tr(Dfip|“) u^n^Pu 

\ OX ox) ox 

u^B^Pu (m^Pu)(u^Dm) 


>0, VccGK'^. 


The positivity of Yp follows from the Cauchy-Schwarz inequality using 
(w^DffP m) 2 = Tr(v^(it 0 u^)^/B VB^VP)^- Note that, for P := D, the 
matrix Yp would simplify to Y from Lemma 12.101 

The differential inequality (12.401) for 5y,(/(t)) then yields the claimed exponential 
decay of 5^, (/(<)). □ 


This concludes the first step of the modified entropy method. In the second step 
we want to prove exponential decay of the relative entropy e^(/(t)|/oo)- In the 
standard entropy method this is achieved by integrating the differential inequality 
(12.311) for Iijj{f{t)) in time, since e'it) = —I^{f{t)). But here, this is not possible, 
since 5'^(/(t)) is not the time derivative of e(t). Instead, we shall use convex Sobolev 
inequalities (cf. §3 of [7]; [39]), which give a simple relation between these two 
functionals. In fact, the functional 5^(/) controls the relative entropy ei/,(/|/oo): 

Lemma 2.13. Let P be some fixed positive definite matrix. Then, the following 
convex Sobolev inequality holds V5 G L)|_(R'^) with g dx = 1: 

(2-41) e^(5l/oo) < ^ 

where both sides may be infinite. The constant Xp > 0 is the smallest eigenvalue of 
P, i.e. 

(2.42) P>ApI>0. 

Proof. As an auxiliary problem we consider the following symmetric non-degenerate 
FPE for g = g(t,x) on L'^{f-^): 

(2.43) dtg = div (/^PV^) , 

Joo ^ 

d I 3 ^ I ^ 

with /oo = (27r)“2e 2“. This is motivated by the fact that S,p{g) is the true 

Fisher information for the evolution under (12.431) . Obviously, we have goo = foe- 
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We also note that (j2.42ll is the (standard) Bakry-Emery condition for p.43|l . since 
its steady state potential is A{x) = \x\’^/2 (cf. p.29|l l. 

Hence, the entropy method implies exponential decay of g{t) to goo with rate 
2Ap (cf. (|2.351) 1. Moreover, the entropy inequality (I2.34|) is already the claimed 
result. □ 

Combining this lemma with Proposition 12 .1 2l readily yields exponential decay of 
the relative entropy, provided that ^^(/o) < oo: 

Theorem 2.14. Assume condition (A). Let ^ generate an admissible entropy, let 
f be the solution to \2.1]) with an initial state satisfying 5'^(/o) < oo, and let fi := 
min{Re{A}|A is an eigenvalue ofC}. Let {Am|l < rn < mo} be the eigenvalues of 
C with Re{Am} = g, and let P be defined as in Lemma \2.11i 
(i) If all \m, 1 < w < mo, are non-defective, then 

(2.44) e^(/(t)|/oo) < ^5^(/o)e-"^‘, < > 0. 

(a) If Xm is defective for at least one m G {!,..., mo}, then 

(2.45) e^(/(t)|/oo) < t > 0, 

for any e G Here, S,p{f,e) denotes the modified entropy dissipation 

functional 1)2.361) with the matrix P = P(e). 

We remark that the multiplicative constant in (12.451) is e-dependent, with S'^(/o, e) 
oo as e \ 0. In (12.441) the exponential decay rate is indeed sharp (cf. §6 of [S]). 
Also, it is independent of the normalizing transformation in (12.31 since the drift 
matrices C and C are similar. But compared to the standard entropy method, the 
above result is not yet fully satisfactory: In the decay estimate (12.351) the initial 
condition is only required to have finite relative entropy. By contrast. Theorem 
12.141 requires the initial state to satisfy 5'^(/o) < oo, and this functional is closely 
related to a weighted iJ^-norm. This “deficiency” of Theorem 12.141 can be lifted by 
exploiting the hypoelliptic regularization of m, cf. also Proposition 12.61 a'). The 
following result is a generalization of Theorems A.12, A.15 in [41] (expressed for 
quadratic and logarithmic entropies) to all admissible i/^-entropies. For its proof we 
refer to Theorem 4.8 in [5]. 

Lemma 2.15. Let condition (A) hold, fo G L^{W^) with fj^dfo dx = 1 and 
e^(/o|/oo) < oo. Let f(t) be the solution of i2. II) with initial condition fo, and let 
T be the minimal constant such that (12.211) (or, equivalently, (I2.22p ) holds. Then 
there exists a positive constant Cr > 0 such that 

(2.46) S^ifit)) < c,t-(2”+i)e^(/o|/oo), Vt e (0,1]. 

With this ingredient we are ready to state our final result: 

Theorem 2.16. Assume condition (A). Let if generate an admissible relative en¬ 
tropy and let f be the solution to with initial state fo G L)(.(R'^) such that 

eV'(/ol/oo) < oo. Let /i := min{Re{A}|A is an eigenvalue ofC}. Let {Am|l <m< 
Too} be the eigenvalues of C with yt = Re{Am}) and let 

e{t) := ev,(/(t)|/oo). 


Then: 
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(i) If all Am, 1 < m < mo, are non-defective, then there is a constant c > 1 
such that 


(2.47) e(t) < ce ^'"*e^(/o|/oo), Vt>0. 

(ii) If Am is defective for at least one m € {1,... ,mo}, then, for all e € 
there is Cg > 1 such that 

(2.48) e(t) < Cee“^(^"®)*e^(/o|/oo), Vt>0. 

Proof. Let P be defined as in Lemma [2.111 Fix some (5 > 0, and let k := /i in 
case (i), and k := — e in case (ii). Using the convex Sobolev inequality (12.411) . 

Proposition 12.121 and Lemma [2.151 we compute loi t> S: 

eit) < ^5^(/(0) < 

(2.49) < 

For t < 6, the monotonicity of e{t) (cf. (12.141) 1 implies 

(2.50) e(<)<e(0). 

Writing cs := e^'^^maxjl, 2 Xp%,-+i I combining (12.491) . (12.501) yields 

e(t) < C5e(0)e“^'‘*, Vt>0. 


□ 

We remark that the exponential decay rate 2k is sharp here, but the multiplica¬ 
tive constant c will in general not be sharp. 

To close this section we shall now briefly illustrate the mechanism of the presented 
modified entropy method. To this end we return to Example 1 2.2 1 and the “distorted” 
vector norm 

\x\p := x/{x,Vx) , 

with P > 0, that was already used in Figure [1] The drift characteristics xft) 
corresponding to (12.191) satisfy xt = —Cx. For the decay of this P-norm along a 
characteristic we obtain 

(2.51) \x\l = -2x'^PCx = -x'^{PC + C^P)a; < -2fj.\x\l, 

where we used in the last step the matrix estimate (12.371) for Q := C and the 
notation /i := min{Re{A}|A is an eigenvalue of C}. So, /r is the spectral gap of 
C, i.e. the distance of cr(C) from the imaginary axis, and it determines the best 
possible decay of x{t). Due to (12.511) . \x\p realizes this optimal decay uniformly in 
time. 

The matrix P determining this “distorted” vector norm is defined via (12.371) , and 
hence it is the same matrix as in the definition of the modified entropy dissipation 
functional S,p{f). □ 
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2.5. Entropy methods for non-degenerate Fokker-Planck equations. 

We remark that the new entropy method from 92.41 is not restricted to degenerate 
FPEs. For non-degenerate FPEs it is in fact a generalization of the standard 
entropy method: Eor a symmetric FPE with constant diffusion and drift matrices, 
the normalization of 92.31 yields D = Cg and D is symmetric positive definite. 
Applying Lemma [2dlji) to Q := C with /r := Aniin(C) admits the choice P := D. 
Hence, Sjp{f) = /^(/|/oo) and the method of 92.41 reduces to the standard entropy 
method. 

For non-symmetric FPEs, however, the standard and modified entropy methods 
differ. Eor regular diffusion matrices D > 0, both methods are applicable and 
yield exponential decay of the solution towards equilibrium. So it is natural to 
compare their performances: Eor applying the standard entropy method to (ED) 
in normalized form (i.e. with A{x) = |a:p/2) we consider the corresponding Bakry- 
Emery condition (12.291) : 

I > AflD-i, 


i.e. A_d > 0 is the smallest eigenvalue of D. Then, §2.4 in (or the analogue of the 
convex Sobolev inequality (12.341) 1 implies exponential decay of the relative entropy: 


(2.52) e^(/(t)|/oo) <e-2^^‘e^(/o|/oo), t>0. 


Note that the multiplicative constant in this estimate is 1. 

For the modified entropy method. Theorem 12 . 141 yields the decay estimate 

(2.53) e^(/(t)|/oo) < ^5^(/o)e-2'^‘ Vt > 0 

in the non-defective case (i), with /j, := min{Re{A}|A is an eigenvalue of C}. For 
the comparison of the two obtained decay rates we have the following result: 


Proposition 2.17. Let the coefficients of a non-degenerate, normalized FPE satisfy 
condition (A). With /i defined above, let {Am|l ^ m < toq} be the eigenvalues of C 
with yt = Re{Am}- Then: 

(i) If all Am, 1 < m < mo, are non-defective, then 

(2.54) 0<AD</r. 

(a) If Am is defective for at least one m € {!,..., mo}, then 

(2.55) 0<XD<ti- 


Proof. For case (ii), let A with Re{A} = ^ be a defective eigenvalue. Let p S 
with IpI = 1 be a corresponding eigenvector, and q G a corresponding generalized 
eigenvector. W.l.o.g. we assume that (q, p) = 0, and q satisfies (AI — C)q = p. 

Next we consider a family of generalized eigenvectors, qs := q-\- Sp, d £ M, which 
also satisfy (AI — C)^^ = p. We compute 

qJ{C -f C'^)qs = qJiXqs-p) + iXqJ -p^)qs = 2Re{A} \qs\^ - 25. 


Using D = Cg and [gap = \q\^ 5-6'^ we obtain for the Rayleigh quotient of D: 

qJ'Dqs 6 


Xd < 


\qs\^ ^ |gP + d2’ 


and (j2.55l) follows for any d > 0. 

For case (i) we only need to replace qs by p in the above computation. 


□ 
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For the non-defective case (i), the inequality (j2.54|) will in general not be strict, 
as can be verified on the following simple example: 

/ 1/5 0 0 \ 

C := 0 1/4 -4 , 

V 0 4 1 y 

with the eigenvalues |, | ± iVl015/8, and D = = diag(|, 1). 

So, the exponential decay rate from the new entropy method is always at least as 
good as the rate from the standard entropy method, but often better. The first rate 
2\d from (12.521) gives an estimate for the local decay rate of the relative entropy. It 
reflects the (in absolute value) smallest slope of the relative entropy at any t > 0. 
More precisely, it is, pointwise in time, a lower bound for the local decay rate, i.e. 

non-symmetric FPEs with linear drift it is well known (cf. §2.4, §3.5 
in 0) that this rate is optimal (as a pointwise estimate). In Figure [S] the initial 
condition is chosen such that the function on the r.h.s. of (12.521) is indeed tangent 
to e{t) at t = 0. 

By contrast, the estimate (12.531) describes the global decay. Hence, its multiplica¬ 
tive constant has to be larger than 1 for non-symmetric FPEs. In some examples, 
the r.h.s. of (12.531) is even the perfect envelope of e{t), see Figure [51 

Example 2.6. We consider the non-degenerate, non-symmetric Fokker-Planck 
equation m with 

D = diag(l/4, 1), C = 

which is normalized. Here we have Xd = \ and /i = |, and the local and global 
decay estimates are shown in Figure [51 □ 

So far, we only discussed the modified entropy method for FPEs with constant 
diffusion and drift matrices. Its generalization to some cases of non-symmetric 
FPEs with non-constant coefficients is the topic of the subsequent chapter. 

3. Kinetic Fokker-Planck equation with non-quadratic 

potentials 

In this chapter we shall illustrate how the modified entropy method from (12.41 
can be extended to kinetic Fokker-Planck equations (1^ with non-quadratic po¬ 
tentials V = V{x) (i.e. a drift term that is nonlinear in the position variable). A 
motivation for the following analysis is its possible application to a future study of 
Fokker-Planck-Poisson equations with a quadratic confinement potential and the 
self-consistent potential acting as a perturbation. Refer to [7l §4.2], for the large 
time analysis of a non-degenerate drift-diffusion Poisson model. 

Several proofs of the entropy- and L^-decay of this equation have already been 
obtained in the last few years: In m. algebraic decay was proved for potentials 
that are asymptotically quadratic (as jccj —>■ oo) and for initial conditions that are 
bounded below and above by Gaussians. The authors used logarithmic Sobolev 
inequalities and entropy methods. In |22| . exponential decay was obtained also for 
faster growing potentials and more general initial conditions. That proof is based 
on hypoellipticity techniques. In m exponential decay in was proved, allowing 
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Figure 5. Entropy decay for the non-degenerate, non-symmetric 
Fokker-Planck equation (I2.1|) with D = diag(l/4, 1), C = [1/4 — 
4 ; 4 1]. Solid red curve: decay of the logarithmic entropy ei{t); 
dotted blue: The estimate of the local decay rate from the standard 
entropy method is tangent at t = 0; dashed black: estimate of the 
global decay rate from the hypocoercive entropy method, (colors 
only online) 


for potentials with linear or super-linear growth. This chapter will now provide 
an alternative proof of exponential entropy decay for (12.2p with a certain class of 
non-quadratic potentials and for all admissible relative entropies e^. 

The kinetic Fokker-Planck equation (I2.2|l has a unique normalized steady state 


(3.1) 


/oo(a:,u) = exp V{x) + 


2 J 


x,v G 


for potentials V(x) with lim|2;|_>oo F(a:) = oo sufficiently fast such that foe G 
^i(]g2n)^ see [40]. An additive normalization constant is included in V. 

We rewrite (12.21) again in the form (12.31) . such that 


( 3 . 2 ) dtf = Lf:=dW^[T>y^f + G{Of], 
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where ^ := {x, v)’^ G d = 2n, D is a block diagonal diffusion matrix and G a 
drift vector field given by 

0 0 


D = 


0 al 


and G{x,v) = 


— V 

VxV + vv 


respectively. 

The positivity of solutions of (12.21) with non-negative initial datum can be proved 
using the sharp maximum principle [23]; see also [5j Proposition 7.1]. 

We introduce the modified entropy dissipation functional S^{f) as in (I2.36L 

Sdf) ■■= [ 

Jr^ ^/oo^ ^ Joo'' ^ /oo'' 


with a positive definite and ^-independent matrix P S to be chosen later. The 
time derivative of S^{f{t)) is estimated as in the proof of Proposition 12.121 — apart 
from not normalizing the equation—and it satisfies 

H f rP" 7^ fP' F' 

(3.3) <-J - R)^P + P^(D + R)]«/oo de, 

where u := E{^) := ^[V{x) -f and R = f ^ ^ ^ 

analogy to 

(3.4) 


ixd 


In 


we define the matrix 

d^E 


Qix) := (D - R) 




0 

a'^v I 




I 

ul 


If we can find an x-independent, symmetric, positive definite matrix P > 0 and a 
constant k > 0 , such that 

(3.5) Q(x)P-bPQ^(x)-2KP>0 VxGM", 
then the right-hand-side of p.3l) can be estimated as 

(3.6) < -2k J u^Pu foo d^ = -2KS^{f(t)). 

R<i 

If additionally k > 0 , this would imply exponential decay of S'^(/(t)). 


3.1. Potential V{x) -with bounded second order derivatives. In this sec¬ 
tion we prove in Theorem 13.91 the exponential convergence of solutions of (12.21) to 
the steady state via the modified entropy method. 

To keep the presentation simple, we shall consider from now on only the ID case, 
i.e. X, u G R (d = 2). Furthermore, we shall consider non-quadratic potentials V{x) 
with bounded second order derivatives satisfying 

(3.7) dyi <72 such that 71 < V'\x) <72 Vx G R. 

To apply the modified entropy method, we need to find a symmetric, positive 
definite matrix P and k > 0 such that (1531) is satisfied. We define 

:= such that Q(x) = 

Then, for potentials V satisfying (13.71) with 71 = inf 2 ;gR V''{x) and 72 = sup^.^^ V"{x), 
condition (13.51) is equivalent to 

(3.8) Q-yP + PQ^ - 2 kP > 0 V 7 G [ 71 , 72 ]. 
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Next we collect the conditions on k and on the coefficients of the matrix P: A 
symmetric matrix P G is positive definite iff its first diagonal element and its 
determinant are positive. Condition ( 13 . 51 ) is linear in P, therefore, we consider— 
without loss of generality—matrices 

( 3 . 9 ) P = f ^ G with det(P) = P22 — P12 > 0 . 

\Pi2 P22) 


For given 0 < 12 and 71 <72, we want to determine k > 0 and symmetric, positive 
definite matrices P such that ( 13.811 holds. The matrix 


Q.yP + PQ^ — 2 kP = 


/ 2 (pi2 - k ) 

V “7 + 2k)pi2 + P22 


-7 + (:/ - 2k)pi2 + P22\ 
2 {-IP12 + k)p22) ) 


is again real symmetric. Hence it is positive semi-definite iff its diagonal elements 
and its determinant are non-negative, i.e. pi2 — k > 0, —7P12 + (v — k)p22 > 0, and 


( 3 . 10 ) 0 < 6(k, 7) := det(Q.^P -f PQ^ - 2 kP) 

= 4 (pi2 - K)i-JPl2 + - k)p22) - (-7 - 2k)p12 + P22f' 

for all 7 G [71,72]. We summarize the conditions on the parameters (pi2,P22,«): 

(Cl) det(P) = P22 - p ?2 > 0 P22 > P ?2 > 0, 

(C2) K > 0, 

(C3) pi2 > k(> 0), 

(C4) <5(k,7) > 0 for all 7 G [71,72], 

(C5) —7Pi 2 + (i^ — k)p22 > 0 for all 7 G [71,72]- 

Remark 3.1. Condition (C3) and a strict inequality in Cond ition |(C4) imply 
Condition (C5) Let, for some fixed (pi2,P22,k), the Conditions (C3)f (C4) hold 
for a 7-interval with interior F. Then (C4) holds on F with strict inequality; 
hence also (C5) holds on F. By continuity (C4) - (C5) then also hold on F. Thus 
(except for the case of F being the empty set) Condition (C5) follows from Condi¬ 
tions 


(C3) (C4) 


1 )}" X R’*' is admissible, if there exist kq > 0 
such that (Cl) - (C5) hold with k = kq and 71 = 72 = 7o- 

, then {pi2,P22) 


Definition 3.2. A pair {pi2,P22) G 
and 7o G 


Lemma 3.3. If {pi2,P22) is admissible for some kq > 0 and 70 G 
is admissible also for all k G [ 0 , kq] and given 79. 


Proof. The Conditions (Cl) - (C3) continue to hold for all k G [ 0 , kq] and given 79. 
The admissible parameters {pi2,P22) define a symmetric positive definite matrix P 
satisfying Q^pP -I- PQ^p > 2K9P. Due to P > 0 , Q^pP -t- PQ^p > 2K9P > 2 kP for 
all K G [0,K9]. Hence also Condition (C4) is satisfied for all n G [ 0 , K9] and given 
79. Since P22 > 0 , Condition (C5) carries over to k G [ 0 , K9]. □ 

We rewrite (5 (k, 7) with respect to powers of 7 as 

( 3 . 11 ) ( 5 (k, 7) = -7^ - (4pi2 - ‘^vpi2 - 2^22)7 - c(k) 
with 

( 3 . 12 ) c(k) := 4 k(jz - k){p 22 - P^) + i’^Pi 2 - P22Y = 4 k(i^ - K)ai + 02 , 


with ai := {p22 — P12) > 0 due to Condition (Cl) and 02 := (^'Pi2 — ^22)^ > 0 . 
The function c{k) satisfies c( 0 ) = c{iy) = 02 > 0 , hence, c(k) is non-negative for all 
K G [ 0 , J2] and monotonically increasing for all k G [ 0 , !■]. 
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Lemma 3.4. Admissible pairs (^12,^22) exist only for k G [ 0 , ^]. 


Proof. Assume (^12,^22) is admissible for some kq > Then 7 can be increased 
until ( 5 (ko, 7 o) = 0 . Due to Lemma [ 3 ^ 0 > 70)) Moreover, 


(9 dr’ 

=-^M=8ai(K0-i) 

and ai > 0 imply kq — < 0, contradicting our initial assumption 


□ 


Remark 3.5. 6 { k , 7) describes a parabola (as the function (13.lip of 7) and i5(k, 7)17=0 
—c{k) < 0 for HI G [ 0 ,u]. Therefore, each 7-interval with 6{k,^) > 0 is either in¬ 
cluded in or in But, in the latter case, V"{x) < 0 for all x G K, which 
would not give an integrable steady state. Hence, only 7 > 0 is relevant. 


Next we establish an important condition: < v. 

Proposition 3.6. Let 0 < 71 < 72 he given. If and only if they satisfy the con¬ 
dition ypfi — ypf l < V, then there exists an admissible pair {pi2,P22) satisfying 
Conditions (Cl) (C5) for some kq > 0 and for all 7 G [71,72]. 

The proof is deferred to Section 13.21 


Theorem 3.7. Suppose 0 < u and 0 < 71 < 72 satisfy ypff — < v. Then 

the following {pi2,P22) G R)]" x R+ are all admissible pairs for K^ax G [ 0 , the 
maximal possible value of k, and for all 7 G [ 7 i, 72 ]' 

(Bl) //371 -1-72 < then Umax = 7 - 5 - 471 and 


(Pl2,P22) = (f + - 371 - 72, - 271 -h rUV^J /2 - 371 - 72)) 


with T G [— 1 , 1 ] satisfy the conditions ( Cl ) ( C5) 

^m.n.pT. — 


(B2) //371 -1-72 > then 


72-71 


satisfy the conditions ( Cl ) ( C5) 


2y/2(yf+j2)TiA 


, P12 = 7 and P22 = 


The proof is deferred to Section [?31 


Remark 3.8. The expressions for Kmax, P 12 andp22 are continuous at the interface 
371 -I- 72 = . 

Following Theorem 13.71 we obtain for given u > 0 and 0 < 71 < 72 < (u -|- 
that a symmetric positive definite matrix P and n = Kmax > 0 exist such that (13.511 
holds. Hence, the modified entropy method yields the following theorem. 

Theorem 3.9. Let '0 generate an admissible entropy and let f be the solution to 
the kinetic Fokker-Planck equation i2.2\} with a potential V(x) satisfying (j3.7p and 
an initial state /o satisfying S.ip{fo) < 00 . Under the assumptions of Theorem \3. 7| 
we then have: 

(3.13) ey,(/(t)|/oo) < cS'y,(/o)e“^'^”‘“"*, t>0, 

for some constant c > 0 independent of /o and Kmax given in Theorem \3.'T\ 

Proof. We already noticed that, following Theorem 13.71 we obtain for given u > 0 
and 0 < 7i < 72 < (u -I- y/^Y ^ symmetric positive definite matrix P and k = 
Kmax > 0 such that (13.51) holds. Consequently, inequality (13.6p follows and implies 
the exponential decay of the modified entropy dissipation functional 

(3.14) S^{f{t))<S^{fo)e-^^\ t>Q. 




















tol-^ 
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Figure 6. A visualization of the (71,72) subset such that for a 
given 0 < u there exist parameters (pii,_pi2,P22, k) G x Rq x 
M+ X [0, |] satisfying conditions (C1)-(C5) according to Theo¬ 
rem [S21 


Moreover, due to Lemma 12.131 the convex Sobolev inequality 

ev.(5l/oo) < , Vg e LI{R^) with [ g = 1 


holds, where Xp > 0 is the smallest eigenvalue of P. Thus (13.1311 follows from 

(Cmi) . □ 
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In a previous work [H §7] the authors considered potentials of the form 

(3.15) V{x) = luq —+ V(x) with ||I/"||j;,oo < oo, and wq 7 ^ 0. 

Following the proof of Lemma 12.111 a matrix P, corresponding to the potential 

O 2 

term ujq can be constructed as 


(3.16) 


and 


P := 


— 2 a;Q 


2 CJg 


if 4wg < , 


if 4wg > , 


(3.17) 


2 ko 




if 4a;^ < , 
if 4a;^ > . 


Proposition 3.10 (O Proposition 7.3]). Let 4ajQ 7 ^ and let V from p.l5l) 
satisfy ||l/"||ioo < |_ 4 (^ 21 defined in (13.171) . Then the modified 

entropy dissipation S-,p{f(t)) with the matrix P chosen in (13.161) satisfies 


S^{f{t))<S^{fo)e ' 



l|v"llr,°° V 


for t > 0. 


Theorem 3.11 ((5] Theorem 7.4]). Let if generate an admissible entropy and let f 
be the solution to the kinetic Fokker-Planck equation h2.^) with an initial state /o 
satisfying S',/,(/o) < 00 . Under the assumptions of Provosition IKJU we then have: 

J.O- 

(3.18) e4fit)\f^)<cS^{fo)e ^ ^ ^> 0 , 


for some constant c > 0 independent of /g. 

The defective case 4 wq = is omitted in [3; but it is noted that a matrix 
P = P(£) could easily be constructed from the proof of Lemma [2. Ill fill. 

To compare the decay rates in Theorem l3.9l and Theorem l3.11l we have to relate 
the parameters in Theorem 13.111 with the parameters 71 and 72 in Theorem 13.91 
Moreover, in Theorem 13. Ill the parameter ujq has to be chosen as to optimize the 
decay rate. 


Proposition 3.12. If Q < v and V{x) with 0 < 71 := 

11 ^ ^ ^ 11 

given, then the largest rate k := sup^^j kq ~ 2 4 21 ° 

to Kmax *71 Theorem 


inf P" < supP" =: 72 are 
in Provosition \3. 10\ is equal 


Proof. For given 0 < u and 71 < V"{x) < 72 , we need to decompose V" as 
V" = Wg + V" such as to maximize the function 


\/|l/2 - 4w^| 


= Ko(wo) 


max{|72 - cugl, I 71 -cugH 


k{uJo) = Ko{oJo) 
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with Ko(wo) given in ( 13 . 171 ) . After distinguishing several cases, one obtains that 
K(a;o) = Kmax for 

for < 371 + 72 , 
for 371 + 72 < . 

□ 



In case 71 =72, the admissible potentials in 93 .II are V(x) = 71^ + Cia; + C2 for 
any constants Ci, C2 G R. Consider the limit 71 —>■ 72 in Theorem 13.71 we recover 
in the limit 71 —>■ 72 the decay rate and matrix P from § 7 ] by choosing r = 0 in 


the case KBl) 


3 . 2 . Proof of Proposition [ 376 l 

Lemma 3 . 13 . Let (pi2,P22) be admissible for some kq > 0 and 70 > 0 . Then 
(,Pi2,P22) is also admissible for kq and exactly for 7 £ [71,72] with 

( 3 . 19 ) 71.2 = + vpi2 +P22T \J (-2pi2 + vpi2 + P22Y - c(ko) > 0 . 


Proof. Conditions (Cl) - (C 3 ) and (I(ko, 7o) > 0 hold, since {P121P22) is admissible. 


Consequently, the equation ^(^0,7) = 0 has (one or two) real solutions 71 < 72, 
satisfying 0 < 70 G [71,72] and Condition (C 4 ) holds for all 7 G [71,72]- Due 
to Remark [331 0 < 71 < 72. Moreover, (5 (ko,7) > 0 for 7 G (71,72) and Condi- 

Condition (C 5 ) follows for all 7 G (71,72), and for all 

□ 


tion 
7 G 


holds. Hence, 

[71,72] by continuity, see R.emark l 3.11 


Remark 3 . 14 . Due to ( 13 . 121 ) . Lemma IHITI and Lemma [ 3. 131 the possible 7-interval 
decreases strictly monotonically with k (as expected from Q-yP + PQ^ > 2 ztP). 
For any fixed z^,pi2,P22, the largest possible 7-interval is obtained for k = 0 , i.e. 
with c(0) = a2 = {vpi2 — P22Y ^ 0- 

Proof of Proposition \ 3 .(^ Following Lemma 13.131 and Remark 13.141 we seek the 
largest 7-interval [71,72] (which maximizes 72 — 71 for fixed pi2,P22) and con¬ 
sequently set K = 0 . The expressions for 71 < 72 in (I 3 . 19 |) and zt = 0 yield 

- 2 p ?2 + '^Pi2 + P22 = , 

2\/(-2p?2 + ^-^12 +P22)^ - c(0) = 72 - 7i , 


or, equivalently with oi\ = P22 — P12 and := pi2{v — P12), 

ai+a3 = {p22 - PI2) + {Pi2{y - P12)) = =: P2 , 

ai as = (p22 - Pi2)(pi2(i^ - P12)) = =: / 3 i > 0 . 

Combining the last two equations, we derive 

—as (/?2 — as) + /3i = 0 

which has two real solutions q;3_± = ^ ± ^ — Ajdi = We recall 

<^3 = Pi2{v — P12), which has real solutions pi2 if and only if as < ■^. Due to 
0 < 7i <72, this restriction is equivalent to 

( 3 . 20 ) 


= V32 ± ^ • 
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For 0 < 7i <72, Condition ( 13.2011 with “+” is more restrictive than with 
Therefore, we consider in the sequel 03,- = \ accordance with the 

key assumption in Proposition 13.61 

Condition is equivalent to ai > 0 . Due to 01+03 = q,^ = 

as - = i( y^- we deduce oi = \{y/^ + > 0 since 0 < 71 < 72. 

Condition (C 2 ) is satisfied due to our choice /c = 0 . Next, pi2(u — P12) = 03__ 
has two real solutions 0 < P12,- < Pi2,+ < v with P12,- +pi2,-i- = u; hence, Con¬ 


dition (C 3 ) holds. Due to our construction starting from (I 3 . 19 L Condition (C 4 ) 


holds. Finally, Condition (C 5 ) follows again from Conditions (C 3 ) - (C 4 ) and 
Remark 13.11 □ 


3.3. Proof of Theorem 13.71 

Lemma 3.15. Let (pi2,P22) he admissible for some kq > 0 . Then pi2 < u — kq. 
Proof. By Lemma [ 3. 131 71^2 > 0 . Hence, the discriminant in ( 13 . 191 ) satisfies 
0 < (-2pi2 + vpi2 + P 22 T - i^Pi2 - P22f - 4Koai(u - Ko) 

= 4ai[pi2(u - P 12 ) - Ko(u - Ko)] 

which is equivalent to kq < P12 < u — kq since oi > 0 . □ 


Remark 3.16. The maximal value of k given by Lemma 15 ^ i.e. k = -I, is possible, 
but only for quadratic potentials: It implies pi2 = f, 7 = 7 i = 72 = P22 > x 


to 


(Cl)). 


Lemma 3.17. For 7,71,72 given as in Theorem let Umax denote the maximal 
decay rate and let V denote the set of admissible pairs (^12,^22) S x R+ (w.r.t. 
the whole interval [71,72]^- Then, 

(a) V is convex and compact; and V lies in the interior of the set defined by 
the inequalities \ (Cl )\ an d \(m\ 

(b) V is a finite, possibly one-pointed, line segment with p 12 € [Pi2>Pi2]- 


Proof, (a) The convexity is clear from (j 3 . 8 |l . 


(C 3 ) and Lemma [ 3.151 imply the boundedness of pi2. ( 13 . 101 ) yields an upper 
bound for P22 (by considering the balance of P22 and the linear terms in ^22)- For 
0 < 7i < 72, no points of V can lie on the curve P22 = P12 (cf- (Cl)), since otherwise 
we would obtain: = —(7 +P12 ~ upi2)^, and|(C 4 ) would only be true for 

a single value of 7. Hence, the strict inequality (Cl) also holds for accumulation 
points of V (for (C 2 ) (C 5 ) this is trivial). This implies that the bounded set V is 
closed. Hence, P is compact. 

By the same argument we have for all (pi2,P22) G P- 


( 3 . 21 ) pi 2 > K, 

since otherwise d{K,j) = — (7 — (u— 2 k)k — p22)^■ Hence, P = P lies in the interior 


of the set defined by the inequalities (Cl) and (C 3 ) 

[(b)] For each fixed (pi2,P22) G P, we have 

(3.22) ^(^max, 7l) — 0 OF ^(^max, 72 ) ~ 0 


(or both): Otherwise, due to Remark 13.141 and ( 13 . 211 ) . Kmax could be increased 
slightly, which contradicts maximality of Umax- For fixed pi2, assume now that 
^Ipi2 ■= {P > 0 I (Pi2,p) G F} is not one point, but rather a closed interval (due to 
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the convexity of V). Then, one of the equations in (|3.22|1 holds for more than two 
values of p 22 - But this is impossible, since 6{k,j) = 0 is a quadratic equation for 
P22 (cf. (I3.8p b Hence, consists only of one point and V is a. line segment. □ 

By Lemma 13.171 V is uniquely determined by its endpoints. 


Lemma 3.18. Let^, ^ 1,^2 be given as in Lemma \3.1'T\ For an endpoint {pi 2 ,P 22 ) ^ 
V we have S{Kmax,ll) = 5{Kmax,l2) = 0 . 


Proof. W.l.o.g. we now assume that 5{Kmax,li) = 0 and 5(Kmax, 72) > 0. So the 
inequalities (C4) for 7 = 72 and (C3) hold strictly, as well as (C5) for 7 = 72 
(due to Remark ITD) . Hence, (Cl) (C5) also hold for 7 = 72 and all {pi 2 ,P 22 ) in 
a small neighborhood of (Pi 2 ,P 22 )- 

Finally we consider, for pi 2 fixed, S{Kmax,Ji) = 0 as a quadratic equation for 
P22- The discriminant for its real solvability reads 

[pi2U - 2k{u - k)+ 71]^ + [-71 + {u - 2«;)pi 2]^ + 4(pi2 - At)7ipi2 . 


For 7 i > 0 this is positive due to (13.2111 . and for 71 = 0 since k < ^. Hence, 
S{Kmax,7i) = 0 is also Solvable for P 22 , if P 12 lies in a small neighborhood of Pi 2 . 
Thus, Pi 2 is not an endpoint of the line segment V. □ 


Proof of Theorem 17.71 Step 1: For given 0 < 71 < 72 , we shall first find admissi¬ 
ble endpoints (^ 12 ,^ 22 ) G V such that (C1) -(C5) hold exactly for all 7 G [ 71 , 72 ] 
with the maximal k € [0, fj. The expressions for 71 < 72 in (13.1911 yield 

- 2 pi 2 + '^Pi2 +P22 = , 


\/ (-2p' 


+^'Pi 2 +^22)^ - c(k) = 


or, equivalently with ai = P 22 — P 12 and 03 := pi 2 {v — P 12 ), 

(3.23) ai+a 3 = {p 22 - PI 2 ) + {pi 2 {v - P 12 )) = =: h , 

(3.24) Qfi [03 - K {v - k)] = =: /3i > 0 . 

For the line = [32 — cti to intersect the hyperbola = ^ + n {v — k) a.X some 
Oi > 0, we require that 0 < k {i' — k) < P 2 , see also Figure [3 
The solutions of (I3.23ll - ()3.24|l read 

(3.25) q;3,± = ± 1 ^^ ^ 

We seek the maximum k G [0, 7 ] such that 03 G K. (for k = 0 this always holds 
by the proof of Proposition 13.611 . This maximal value is either obtained as k = ^ 
or when the discriminant of (I3.25|l is zero. The latter case implies 

-K [v- k) = 2 -132 = -71 • 

This is solvable (for k) in R iff yielding /c = .| — \J^ G [0, ^]. Hence, 

for the solvability of (I3.23II - (I3.24II in R, we obtain 

for y/^ < I, 
for y/^ > I . 
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as 



Figure 7. For the line as = 132 — ai to intersect the hyperbola 
as = ^ + (3s aX some ai > 0 , we require that 0 < (3s < (32- 


Using K in (I3.25|) yields one or two values for as > 0. Next, we need to check the 
solvability of as = Pi 2 (u — P 12 ): To obtain pi 2 € M, we must have 

(3.26) as < 4 ■ 


Since 03 with the negative sign gives the weaker constraint, we shall use only 03 ^- 
in the sequel. Now, we have to distinguish between three cases: 


(AI) 

(A2) 


< I and 371 + 72 < The unique as-id) = = 3ji+j2 

satisfies condition (13.261) . Hence, (I3.23l) - (l3.24p yield the two endpoints for 


(pi 2 ,P 22 ) given in (Bl) 


^ and 371+72 > v^. Here as-iji) violates condition (13.261) . Hence, 

i^max has to be chosen smaller than k. Since as^-{K) is monotonically 

2 

increasing, the obvious choice as ■= ^ also yields the maximal value of n: 
Equations (I3.23I) - (I3.24I) give ai = (32— 

[ 0 , - 1 ] in case [(B 2 ) 


and hence K (u-k) = 


with the solution 


+ (3s 
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(A3) Using f yields from p.25l) 

«3.- = ^ + 4 - W(/32-^)2-4/3i > 0. 

But one easily checks that it violates again condition (13.261) . As in case 
(A2) one chooses := and the expressions for (pi 2,^>227 


case 




(B2) follow. 


Finally, Conditions (Cl) (C5) are easily verified for each subcase. 

Step 2 : The whole interval of solutions in (Bl) is obtained due to the convexity 
oiV. □ 


4. Fokker-Planck equations with non-local perturbations 

4.1. Introduction. In this chapter we investigate properties of the following 
class of perturbed Fokker-Planck equations: 

(4.1a) /* = V • (DV/ -b Cx/) -b 0/ = L/ + 0/, 

(4.1b) /(t = 0,x) = (^(x). 

Thereby / = f(t, x), and t > 0 and x G R", with n G N. The matrices D~^C, D G 
are symmetric and positive definite, hence L is a symmetric Fokker-Planck 
operator in L^(K";exp(ix^D~^Cx)), in fact it is a special case of (12.41) . The 
perturbation is given by a convolution 0/ = 'd*/ with respect to x. The convolution 
kernel i? is assumed to be t-independent, and massless, i.e. i?(x) dx = 0. To 
keep the solution / real valued we shall consider here only real valued kernels "d, but 
the analysis would be equally valid for complex i9’s. Further, technical assumptions 
are specified in the beginning of Section l4^ 

The aim of this chapter is to make a spectral analysis of the perturbed Fokker- 
Planck operator in an appropriate weighted L^-space, and to show the existence of 
a unique (up to normalization) stationary solution. Furthermore, the exponential 
decay of any solution of 63) to the stationary solution is proven. 

The following analysis is structured as follows. After notational preliminaries in 
Section IQ we investigate in Section 14.31 the unperturbed Fokker-Planck operator 
in several functional spaces. First, we recall some of its properties in the L^- 
space weighted with the reciprocal of the zero eigenfunction (this weight grows 
super-exponentially), in which the Fokker-Planck operator is self-adjoint. Then, a 
spectral analysis in a larger, exponentially weighted space is given for this operator. 
Finally, in Section lT^ we consider the influence of the pertnrbation 0 on the spectral 
properties of the unperturbed Fokker-Planck operator in the exponentially weighted 
space. 

Equation (gH]) is a toy model for the Wigner-Fokker-Planck equation, see [5]. 
Other examples for equations of this form can be found in gD] and [53]. The 
following analysis of (14.11) is a generalization of the results published in |37j , where 
only the case C = D = I was considered. In this chapter we use a similar approach 
for proving the desired results. However, several proofs and technicalities differ 
from [37] . 

4.2. Preliminaries. We use the convention N = {1, 2,...}, and we write No := 
N U {0}. Given a complex number z € C the complex conjugate is denoted by z. 
For n G N the elements of C” are denoted by bold lowercase letters. Given some 
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vector z S C", the z-th component is denoted by Zi, and we write z — [zi,..., Zn]'^ 
as a column vector. For a multiindex k £ Ng we use the notation z^ := z^^^ ■ ■ ■ 
Given a real number s > 0 we dehne 

For i € {1,..., n} the z-th unit vector in C" is denoted by e^. For every 1 < p < oo 
we define the corresponding p-norm on C" by 

n 1 

|z|p := (X! l<P<oo, 

i^l 

|z|oo := max \zi\. 

l<i<n 

With respect to the norm | • |p the open ball in C" with radius r > 0 and center 
a G C" is defined by 

BP {a.) = {z G C" : |z — a|p < r}. 

Its complement in C" is denoted by := C”\i?^(a). Whenever we work in R” 

instead of C" we use the same notation. Matrices are denoted by bold capital letters. 
For a matrix M G and a real number s > 0 we define := exp(Mlns), 

using the matrix exponential. 


On a domain C R” we call a real-valued function w G a weight 

function if G The corresponding weighted L^-space Lf{VL-,w) is the set 

of all measurable functions such that the norm 

\\f\\<:i;w ■■= |/(x)pw;(x)dx) " 

is finite, and the corresponding inner product is denoted by (•, ■)q-^w 

Also, we introduce weighted Sobolev spaces. For two weight functions wq and wi 
the space wo,Wi) consists of all functions / G wq) whose distributional 

first order derivatives satisfy dffdxj G L‘^{V,-,wi) for all 1 < j < rz. We equip the 
space H^(fl;wo,Wi) with the norm 



which makes it a Hilbert space, see Theorem 1.11 in [55]. If H = R" we shall omit 
the symbol H in these notations. We call two sets of weight functions equivalent 
if the corresponding weighted spaces are the same. In the case where the weight 
functions are equivalent to the constant function, we omit the weight function in 
the notation, e.g. L^(H; 1) = 


For functions / G L^(R") we define the Fourier transform of / as 

= [ /(x)e--«dx. 


We use the same notation for the natural extension of the Fourier transform to 
tempered distributions / G o5^'(R"). With this scaling we may identify /(O) with 
the mass (or mean) of /. For a tempered distribution / G J^'(R") and a multiindex 
k G Ng we define 


VV(x) 


dMif 


dxl^ ■■■dxl. 


as a distributional derivative. 
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Furthermore, we present some definitions and properties concerning linear op¬ 
erators and their spectrum. Let X, X be Hilbert spaces. If X is continuously and 
densely embedded in X we write X ^ X^ and X X indicates that the em¬ 
bedding is compact. Given a subset F C X, the closure of F in X is denoted by 
either F or c\x F. “^(X) denotes the set of all closed operators H in X with dense 
domain D{A). The set of all bounded operators A: X ^ X is ^{X, X)\ if X = X 
we just write 3S{X). Thereby || • ||.^(x) denotes the operator norm. For an operator 
A G “^(X) its range is ranH and its null space is kerH. Note that there always 
holds kerH C D{A). A closed, linear subspace F C X is said to be invariant 
under A G ^{X) (or A-invariant) iff D[A) flF is dense in F and ranAly C F, see 
e.g. [I]. For any C G C lying in the resolvent set p{A), we denote the resolvent by 
Ra {0 '■= (C ~ A)~^. The complement of p[A) is the spectrum cr(A), and (Jp{A) is 
the point spectrum. For an isolated subset a' C o'(A) the corresponding spectral 
projection Pa, a' is defined via the line integral 

(4.2) PA,a' :=^^£RA{C)dC, 

where F is a closed Jordan curve with counter-clockwise orientation, strictly sepa¬ 
rating a' from a{A) \ a', with a' in the inside of F and <j{A) \ a' on the outside. 
The following results can be found in Section III.6.4] and [3H1 Section V.9]: 
The spectral projection is a bounded projection operator, decomposing X into two 
A-invariant subspaces, namely ranP ,4 g.' and kerP 4 This property is referred to 
as the reduction of A by Pa, a-’ ■ A remarkable property of this decomposition is 
the fact that cr(A|ranP^ ,,,) = cr' and cr(A|kerP^ ) = cr(A)\CT'. Most of the time 
we will be concerned with the situation where cr' = {A} is an isolated point of the 
spectrum. 

A final remark concerns constants occurring in estimates: Throughout this chap¬ 
ter, C denotes some positive constant, not necessarily always the same. Dependence 
on certain parameters will be indicated in brackets, e.g. C(t) for dependence on t. 

4.3. Analysis of the Fokker-Planck operator. In this section we investigate 
the (unperturbed) Fokker-Planck equation 

(4.3) ft = X- (DV/ + Cx/). 

Indeed we can find coordinates that simplify this equation. To this end we pro¬ 
ceed similarly to the “normalization” of the Fokker-Planck operator after Theorem 
12.81 Since D is symmetric and positive definite we may introduce the coordinate 
transformation y = Vb- X. With g{y) := /(x) equation (14.31) transforms to 

(4.4) = Vy . (Vyff-H Cyg), 

with c = Vb cVn. Since C is symmetric and positive definite, we may express 
the variable y in terms of an eigenfunction basis of C. Applying this change of 
coordinates to (14.41) yields an equation of the same form, but now the matrix C is 
diagonal (compare to the situation in (12.281) 1. 

Therefore, without loss of generality we shall always assume that D = I, and 
C is diagonal in the following, i.e. C = diag(ci,..., c„) with the entries 0 < ci < 
C 2 < • • • < c„. We introduce c := [ci,..., c„]^. The unperturbed Fokker-Planck 
operator L is then 


L = A + x^CX + TvC. 
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Note that the perturbation 0 in (14.11) still is a convolution in the new coordinates. 
One can check that 

^ := exp(—ix^Cx) 

is a steady state of (14.31) . i.e. a zero eigenfunction of L. The natural (self-adjoint) 
setting for L is the space H := L^(l//i), with the inner product denoted by (•, ■)h- 
There, L is properly defined as the closure of T|c^(b")- This procedure also yields 
the domain D{L). The behavior of L in i? is well studied (cf. [31] [TTl [211135] ), we 
list its main properties in the following theorem. For the case C = I an analogous 
result has been published in m- A complete proof of the following theorem can 
be found in [36] . 

Theorem 4.1. The Fokker-Planck operator L in H has the following properties: 
(i) The operator L = c\h L\c°° on the domain D{L) is self-adjoint and has a 
eompaet resolvent. 

(ii) The speetrum consists entirely of isolated eigenvalues and it is given by 
CT(L) = {-c-k;keN((}. 

(Hi) The zero eigenspace is spanned by /io(x) := det(C/(27r))^/^ exp(—ix^Cx), 
and for every k G Ng the function /rk(x) V^p,o{x) is an eigenfunction 
to the eigenvalue —c • k. 

(iv) For every ( G <j{L) we have ker(C — L) = span{/Xk : C = • k}. 

(v) The family of eigenfunctions {/ik : k G Ng} is an orthogonal basis of FI. 

(vi) L generates a Co-semigroup of contractions (e‘^)t>o in FI, and 

l|e‘^|HjU(H) =e-'==^‘,fcGNo, 

where ci is the smallest entry of c, and FFk '■= span{/rk : |k|i < k — I}"*". 
The following result is useful in the subsequent analysis: 

Lemma 4.2. For every k G Ng the eigenfunction fj,]^ is of the form 

n 

(4.5) Aik(x) = ^o(x) J|p^"(aij), 

where p^^ {xj) is a polynomial of order kj. 

Proof. We prove this by induction. For k = 0 the statement clearly holds true. Let 
it now hold true for some k G Ng, and we deduce the validity for k -|- e^ for any 
^ G {1,..., n}. According to the property /rk = and the induction hypothesis 

we have 

n 

Mk-fe, (x) = de (/ro (x) pf {xj )) 

1=1 

= (/io(x) ( - Cf,xipf{xi) Fpf{xf)'). 

We define the new polynomial p\^^^{xi) := —cixip^'^ [xi) + {p\^{xt))' and it is 
obviously of order -|- 1, since q > 0. This proves (ITB . □ 

For the subspaces Flk, k G Ng, which were introduced in Theorem 14.11 (f^ we 
find the following characterization: 
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Lemma 4.3. Let k € No- There holds / G Hk iff 

(4.6) f /(x)x*^dx = 0, V|k|i</c—1. 




Proof. For this we will rely on the representation ()4.5I1 for the 
then shown by induction. Clearly, we have Hq = H and for k = 


Hi= Ho = 


{feH: [ /(x)dx = 0|. 


Let us assume now that gH) holds for some fc G Nq. According 


Hk+i — 


P n 

{/e-fffe: / /(x) ]^p^" (a::j)dx = 0, V|k|] 


/ik- The result is 
1 we obtain 


to (14.51) we have 

=4 


For f € Hk and |k|i = k we get due to the induction hypothesis 


0 = 


i=i 


/(x)x‘" dx, 


where Ok ^ 0 is the leading coefficient of the polynomial in the integral. All other 
parts of the first integral vanish due to the induction assumption (14.61) . Since this 
holds for all |k|i = k this proves the desired condition for / G Hk+i- □ 


For every k G Nq we define the projection operator 11^ k corresponding to /rk by 
the orthogonal projection 

nL,k := II ■ 

WMh 

With this, the spectral projection corresponding to an eigenvalue C = —c-k is given 
by the orthogonal sum 

nL,c := ^ n^.k- 

keN" 

-c-k=C 

So far we have discussed the operator L in H. However, for investigating the 
perturbed Fokker-Planck operator L + 0 the space H is not convenient. This can 
be illustrated in the one-dimensional case with 0/ := /(x -I- a) — /(x — a), for any 
a > 0. There one can explicitly show that the zero eigenfunction of L -f 0 does 
not lie in H, for more details see [36]. Thus we are forced to investigate L -f 0 in a 
weighted L^-space with a weight which grows more slowly than 1/ho as |x|i —>■ oo. 
It turns out that 

n 

(4.7) c.(x) := ^ cosh jdxi 

i=l 

is a convenient weight function. Thereby /3 > 0 is an arbitrary constant which is 
not yet specified. Note that this differs slightly from the weight function chosen 
in m- However, this choice is more practical for the subsequent analysis. In the 
following we analyze L -f 0 in the weighted space Tl := Lf(u}). The natural norm 
and the inner product in TL are denoted by || • and (•, ■)ui, respectively. 

The space TL possesses a useful characterization via the Fourier transform. 
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Proposition 4.4. There holds f € TL ijf its Fourier transform / possesses an an¬ 
alytic continuation (still denoted by f) to the open set ri ^/2 := {z € C” : | Imz|i < 
/3/2}, with the property 

(4.8) sup |l/(-+ib)||i 2 (R„) < oo. 

bGR" 

|b|i</3/2 

In this case we have: 


(i) For every b G R” with |b|i < /3/2 there holds 

(4.9) /(^+ ib) = J'[/(x)exp(b-x)](^), ^ G K”. 

(a) For every b G K” with |b|i = /3/2 we define /(^ + ib) := J^[/(x)exp(b • 
x)](£), which lies in L'^fEF). With this there holds b i—>■ /(• + ib) G 

See Theorem IX. 13 in [34] for a very similar result. For a detailed proof see 
[36]. Often we shall use the following norm, which is equivalent to || • ||a; due to the 
Plancherel theorem: 


(4.10) 


\\\f\\\l :=E||/( 




2 

L2(Rn) 



2 

L2(Rn)’ 


A useful property of H is the validity of the following Poincare inequality: 


Lemma 4.5. There exists a constant Cp > 0 such that for every f G H^{uj,uj) 
there holds 


(4.11) ll/L<Op||v/L. 

Proof. For this we use the norm ||| • 1^,,. We compute 


iiiwiii5 = EE 


+ i 




|2 


+ 11(0 


L2(R") 


n 

- E (iKo + i2)/(^ + ^2®^)! 


1=1 


2 

L2{E^ 


+ ll(0-if)/(4-i|e£)||?,2(Rn)) 

n 

- (i) E (llA^+ ^2®^)IL2(Rn) + ||/(^ - 


L2(Rn 


= (§) 


= 1 

2 ,„ .,,,2 

LJ ’ 


This proves the Poincare inequality with the constant Cp = 


□ 


Using the above properties of H we can investigate LinH. The following theorem 
is the main result of this section and describes the (unperturbed) Fokker-Planck 
operator in H: 

Theorem 4.6. Let a;(x) be the weight function defined in gll) for any fd > 0, and 
FL := L^{u}) is the corresponding weighted space. Then the Fokker-Planck operator 
A|cj°(R") closable in H, we write C := cln L\(y^(^ny In FL the operator C has 
the following properties: 
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(i) The resolvent of C is compact, and a{C) consists entirely of isolated eigen¬ 
values. 

(ii) The spectrum of C is given by 

a(£) = {-c-k:kGN”}, 

where c is the column vector containing the diagonal entries ofC. 

(Hi) For every A G o'{£) the eorresponding eigenspaee of C is given by 

span{/^k ; k G Nq A —c • k = A}, 

where the eigenfunctions were introduced in Theorem \4. 1\ 

(iv) For every fc G Nq the following is a closed subspace ofTL: 

'^k := |/ G "H : J /(x)x'^ dx = 0, Vk G Nq with |k|i < k — l|. 

TLk is C-invariant, and cr(£|^^) = {—c • k : |k|i > k}. There holds the 
identity 

TL ='Hk® span{/Xk : |k|i < fc — 1}. 

(v) C generates a Co-semigroup of bounded operators {e*^)t>o on TL. For every 
fc G No there exists a eonstant Ck sueh that 

(4.12) ||e‘^|wJU(«) < Vt > 0. 

The rest of this section is dedicated to proving Theorem 14.61 The proof is 
structured into several lemmata and propositions. To this end we begin by showing 
that the Fokker-Planck operator can be defined as a closed operator in TL and we 
characterize its domain. The first preparatory result is the following lemma, which 
is also essential for showing the compactness of the resolvent of the Fokker-Planck 
operator in TL. 

Lemma 4.7. Let Re C > 4(1 -1-/3^-l-Tr C), and f,gG C§°(M.'^) such that {( — L)f = 
g. Then there exists a constant C > 0, independent of f,g, such that 

(4.13) ||/IU+||V/L<C||5L. 

Thereby w{x) := (1 -|- |x| 2 )a;(x). 

Proof. For f = 0, g = 0 (14.131) holds trivially. For / ^ 0 we apply {■,f)t.j to 
(C — L)f = g, and compute 

Re/" gfujdx = Re f (C/- V • (V/-f Cx/))/a; dx 

= ReC / |/|^a;dxH-Re j (V/+ Cx/) - (a;V/+/Vo;) dx 
= l|V/||5 + i [ |/|^(2ReCa;-Aa;-a;TrC-tx^CVw)dx 

2 jR-n 

(4.14) =||v/||2+i/' |/|2^dx. 

^ JR^ 

Thereby we temporarily define z/(x) := 2ReCa; — Aw — wTrC -|- x^CVw. We 
observe that Aw = /3^w and x^CVw = (3 CiXi sinh/3xi > 0 for all x G K". So 
if Re C > ^(1 -I- /3^ -I- Tr C), the function v(x.) is a weight function with z/(x) > w(x) 
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on K". Next we apply the Cauchy-Schwarz inequality to the left hand side of (|4.14|1 . 
which yields 

l|V/|l^ + i||/||^<||/L||gL. 

We now use the Poincare inequality on the first term and u(x) > w(x) on the second 
term, and divide by ||/||(^: 

I|V/L + ||/IU<C||5L. 

Finally we observe that for any fixed ReC > ^(1 -I- /3^ -f Tr C) there is a constant 
C > 0 such that u(x) > Cw(x) for all x G R”. This concludes the proof. □ 

Before we properly define the Fokker-Planck operator as a closed operator in TL, 
we need the lemma below. It determines all formal eigenfunction of the Fokker- 
Planck operator, i.e. the eigenfunctions of the distributional Fokker-Planck operator 
£ in T-L. Thereby, we define the distributional Fokker-Planck operator as £ := 
A-tx^CV-f TrC in the sense of tempered distributions. £ is then a well-defined 
linear map from H into defined on the whole space TL. As a consequence of 
the following lemma it will be straightforward to determine the spectrum of the 
Fokker-Planck operator in TL. 

Lemma 4.8. The distributional Fokker-Planck operator £ satisfies the eigenvalue 
equation £/ = (f for some ^ G C and some f € TL \ {0} iff C ^ • k ; k G Nq }. 

For such values of C, there holds f G span{/Xk : — c • k = C}. 

Proof. Since all the functions are eigenfunctions of L and lie in TL it is clear 
that they are also eigenfunctions of £. In order to show that they already span all 
eigenspaces we consider the Fourier transform of {fi — £)/ = 0 for any ( G C, which 
reads 

(4.15) (C + ieil)f + ^^CVf = 0. 

Now we are looking for / G 'H and ( € C satisfying this (eigenvalue) equation. This 
means that we are interested in solutions / which are analytic in n^/ 2 . Expecting 
/ to be generated from fio by repeated differentiation (see Theorem 14.11 diul) 1. we 
make the ansatz / = p/lo, with p analytic in Tlp/ 2 - This is admissible (and not 
restrictive) since {1q is nonzero and analytic in ^pj 2 . We know that fio satisfies the 
zero eigenvalue equation +€^CV/lo = 0, so after inserting / = pfio in (14.151) 

we obtain the following equation for p: 

(4.16) = -Cp. 

To solve this first order PDE we consider its characteristics: We introduce the 
(unique) solution ^{t) of the ordinary differential equation ^ with ^(0) = 

G C”. It is verified by application of the chain rule that for any such curve and 
any differentiable function p we have 

^p(i(0) = m'^cvpm)- 

In particular, any solution of (I4.16|) fulfills the ordinary differential equation 

-^p(^(i)) = -Cpi^W) 
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along these curves, and it follows Using the fact that |(t) = 

and introducing s = e* (with s S 1 R+) we obtain |(t) = 5*^40 (see Section 1121 
concerning the notation), and consequently we obtain 

(4.17) p(s‘^lo) =P(^o)s“‘^- 

Now p needs to be analytic in n^/ 2 - So (14.1711 implies that Re^ < 0 is necessary, 
otherwise p would have a singularity at the origin ^ = 0 (corresponding to s 0 ), 
which is a contradiction. By induction we deduce from (14.1611 that for all k £ Nq 

i^CV{v'^p) = -(C + c • k)V‘^p. 

Since all derivatives need to be analytic in r 2 ,g /2 as well, the above argument 
proves that either Re C < —c • k for all k £ Nq (which is impossible since C > 0 ) or 
= 0 in 17^/2 for some k £ Nq. So p has to be a polynomial, and we make the 
ansatz 

P(l) = XI 

ksN" 

where pk = 0 for almost all k £ Ng. We now insert this in (14.1611 and obtain 

X (‘^ ■ = -c X 

k 6 N" ksN" 

This holds true iff C = — c • k for all k £ Ng for which pk 7 ^ 0. This proves the first 
statement of the lemma. 

From the above analysis we conclude 

= { X Pk€'")Ao(l)- 

kGNj 

c-k= —C 

Now recall from Theorem 14. II (Hull that pk = holds for all k £ Ng. Hence, 

/ £ spanjpk : — c • k = ^}. So we conclude that the eigenspaces of £ in "H are 
precisely spanned by the pk- D 


Now we can properly define the Fokker-Planck operator in the space 'H. 


Lemma 4.9. The operator L\c’^ is closable inTl, and C := cIt^LIc”- The domain 
is D{C) = {/ £ "H : £/ £ TL}, and for f £ D{C) we have Cf = £/. 

The following proof is based on the proof of Lemma 2.6 in m- 


Proof. According to (I4.14p we have that (L — C)lc“ is dissipative in H if Re^ > 
i(l + + TrC). This implies (cf. [321 Theorem 1.4.5 (c)]) that (L — C)|cj° and 

consequently also L|c~ is closable in H. 

Now we define £ := cI-h £|c“- The domain D{C) consists of all / £ "H for which 
there exists some 5 £ "H and a sequence (/n)neNo C C(j”(R") such that 


(4.18) 


lim ||/„ - f\\^ = 0 , 

n—^oo 

lim \\Lfn- g\\u: = 0 . 


This also implies that {{f—L)fn)neno is a Cauchy sequence in "H. Thus, according 
to (14.1311 (V/„)„gNo is a Cauchy sequence in H. So altogether, (/„)„gNo is a Cauchy 
sequence in the Hilbert space oj). But since we already know that /„—>■/ in 

"H, this implies that even / £ H^{uj,uj). Next we temporarily introduce the weight 
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uj 2 (x) := uj(^) and the corresponding weighted space "^2 := Due to the 

previous results (x^CV/„ + TrC/n)„gNo ^ Cauchy sequence in 'H 2 - According 
to (14.181) . {Lfn)nGNo is also a Cauchy sequence in 712■ Altogether, this implies that 
(A/„) is a Cauchy sequence in ^ 2 - Applying the Fourier transform and the 
norm (14.101) we have that, for every € S {1,..., n}, the two sequences 

are Cauchy sequences in But we also know that /„(• ± if e^) converges to 

/(•±if e^) in L^(]R"). Thus it is clear that A/ G 712, and A/„ —>• A/ in 772 and also 
Lfn —>■ £/ in 772- According to (j4.18l) £,f = g m 772, and since g G 77, we conclude 
that Lfn &f in 77. This proves the inclusion D{C) Q {f G 77 : £,f G 77}. 

Finally we prove that this inclusion indeed is an equality. First we note that 
D{L) C D{C) since L = c\h A|c“ and H ^ 77. So we have the inclusion L C £ for 
the graphs. Let us then take ^ > 0 so large that the estimate (j4.13l) holds. As we 
have mentioned in the beginning of the proof the operator (£ —C)!!?” is (uniformly) 
dissipative in 77, and from Theorem 1.4.5 in [53] it follows that the closure, C — 
is also (uniformly) dissipative. In particular it is injective and thus invertible. So 
(C —£)“^ exists. Now according to Theorem l4.1l d — L: D{L) —>■ iL is a bijection, so 
ran(C —£) D H, which is dense in 77. Due to this and the estimate (14.131) (C —£)“^ 
is a densely defined bounded operator in 77. But by definition — C)~^ is already 
closed, so ran(( — C) =77 and C, G p{C) (and thus p{C) ^ 0). 

For the proof by contradiction we take now this C G p(£)) and assume there 
exists some f* G 77 \ D{C) such that /* G 77. Hence also {(, — £)/* G 77. Since 
C G p{C) we have (C — £)“^(C — £)/* G D{C). Since D{C) is a linear space we have 
:= (C — £)“^(C — £)/* — f* G 77 \ D{C) with {C, — £)/* = 0. But according to 
Lemmawe know that C, G p{C) cannot be an eigenvalue of £ in 77. So /•* = 0, 
contradicting g77 \ D{C). Hence we conclude D{C) = {f g 77 : Zf G 77}. □ 

Lemma 4.10. For any f G p{J~.) the resolvent (^ — L)~^ is compact in 77. 

Proof. We fix (( > 0, and first show the result for this given (}. Choosing (} large 
enough we can apply Lemma 14.71 which proves that (C — C)~^ is an element of 
3§{77,H^{w,uj)). Note that this requires the density of (^“(K") in 77, which is 
assured by Lemma IA.2I in the Appendix. 

Now we shall show that {vu, oj) is compactly embedded in 77. By the definition 
of w (in Lemma EJj) it is clear that for all n G No there holds 

w(x) 1 
sup —= ——, 

|x| 2 >n ^(x) l + n 

which tends to zero as n —> oo. Thus we can apply Lemma lA.ll in the appendix, 
which proves the compact embedding H^(tj7,uj) 77. Hence, the resolvent — 

: 77 ^ 77 is compact. Finally we remark that, according to Theorem 111.6.29 
in [26], the compactness of (C — C)~^ follows for all other f G p{C). □ 

Corollary 4.11. The spectrum a{C) consists entirely of eigenvalues, and a{C) = 
{—c-k : k G Nq}. The eigenspace corresponding to the eigenvalue ( G cr(£) is given 
by span{/ik : C = • k}. 

Proof. We apply Theorem 111.6.29 in [26] which states that (t(£) consists entirely of 
eigenvalues, and the corresponding eigenspaces are finite-dimensional. According to 
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LemmallSlthe eigenfunctions of £ in D{C) are precisely the (formal) eigenfunctions 
of £ in T-L. With this, Lemma [4.81 concludes the proof. □ 


We introduce the closed subspaces Hk H for every k € Nq, which we dehne 
as Hk ■= cluHk, where the subspaces Hk were specified in Theorem 14.11 The 
following lemma gives a characterization of the spaces T-Lk, compare Lemma l4.3l for 
an analogous result in H. 


Lemma 4.12. For every A: G No there holds 

(4.19) nk = {f&n:[ /(x)x*"dx = 0, VkG NO mt/i |k|i < fc-l}. 

Proof. We start from the characterization of the Hk in Lemma [4.31 Our plan is to 
apply Lemma IA.5I in the appendix. For every k G N)) we define the functional 

77k: "H-4 C: / [ /(x)x’^dx. 

We first prove the continuity of the 77 k ■ For k € Ng and f GJi we have 


f /(x)x‘"dx < [ |/(x)a;(x) 


< 


1/2 I 
Rn W(X) 


;(x)i/2 


dx 


dx 


Since w grows exponentially in every direction it is clear that the last integral on the 
right hand side is finite for every k G NO. Thus the 77 k are bounded linear functionals 
in FL. Next we shall verify that the family {r/k : k G NO} is linearly independent. 
If the family would be linearly dependent, there would exist a polynomial p(x) ^ 0 
such that 

[ /(x)p(x) dx = 0, V/ G n. 

JR" 

But this implies p = 0, since C0°(R") C FL. 

Now we have verified the assumptions of Lemma IA.5I Since 

Hk= Pi ker77k|_R, 

|k|i<fc-l 

we conclude that 

FLk ■= cluHk = P ker77k. 

|k|i<fc-l 

The intersection on the right is exactly the set (14.191) . □ 


Corollary 4.13. For fc G Nq there holds the identity 
(4.20) nk = {f Gn:V^Ko)=0, V|k|i<A:-l}. 

Proof. This follows immediately from the fact that tor f G FL and k G NO 

[ xV(x) dx = J-[xV(x)](0) = il*^lw‘^/(0). 

jR" 

We use this in (14.191) and the result follows. □ 

At every A G cr(£) the resolvent map ( i-A Rc{C) has an isolated singularity. We 
denote the corresponding spectral projection of £ by n£_A, which satisfies (14.21) . In 
particular there holds II^^a = dn nL,A, as we will see in the following. 
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Proposition 4.14. For every k G No we have the following facts: 

(i) The space TL can he written as the following direct sum: TL — Hfe©span{^k : 
|k|i<fc-l}. 

(ii) Both spaces TLk anrfspan{/ik : |k|i < k—1} are closed inTL and C-invariant. 
In particular a{C\'Hk) = {~c • k ; |k|i > k}. 

Proof. Step 1 (decomposition of Hk): In H there holds for any fixed k G N 

(4.21) = span{^k : |k|i < fc - l}, 

and for every A G cr(L) we have for the corresponding spectral projection 

(4.22a) ranfli^A = span {/Tk : —c • k = A}, 

(4.22b) kerllL^A = spanj/Tk : —c • k ^ a}. 

For a given fc G N we define the set 

CT/c := {-c • k ; |k|i < fc - 1} C % , 

which is the set of all eigenvalues which “contribute” to (note that there may 
be k G Nq such that —c • k G Cfe but |k|i > k). From (I4.22al) we conclude that 

y ran Hz,, A D . 

Taking the orthogonal complement of this relation yields: 

(4.23) f| kerEzaCFIfe. 

Next we investigate which eigenfunctions need to be added to the left hand 
side of (j4.23(l such that the corresponding span equals Hk- First we observe that, 

according to (I4.22L there holds ^k G (nAecrt kerllz^A) iff Mk G ranllz^A for some 
A G CTfe. This is also equivalent to the condition — c-k G cr^. To complement the left 
hand side of (j4.23L we also require /Xk G Hk, which gives the constraint |k|i > k, 
see (I4.21|l . Hence, we conclude that 

(4.24) Hk = (^ f] kerHi^A) ©_l span{/Xk : -c • k G cr^ A |k|i > k}. 

X^CTk 


Step 2 (decomposition of TL): For ( G p{C) we have Rl{0 C Rc{0 (in the sense 
of graphs), and as a consequence the spectral projection for A G cr(£) satisfies 
Hi, A Cl Hzi^A; see (14.2D . Furthermore, both Hz^^a and H/i.a are bounded projections 
in H and R, respectively. Due to Lemma lA. 61 in the appendix there holds 

(4.25) kerHc^A = clzz kerHz.A and ranH^^A = clzz rannz,A. 

Since the projections are bounded we have TL = kern£ a © ranH/; a, and both 
components of the direct sum are closed subspaces of R, see Section III.3.4 in [26] . 

Step 3 (decomposition of Rk): Due to the arguments of Step 2 we obtain, by 
applying the closure in R to (I4.24D : 

(4.26) Rk = (^ y kern£_A^ © spanj^k : —c • k G ct*, A |k|i > k}. 

X^(Tk 
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Notice that au is finite. The sum is still a direct sum, since every in the “span- 
term” of the right hand side lies in the range of some n£^A with X G ak- Altogether 
this implies that Hk is a closed subspace of H such that 

T-L = 'Hk® span{^k : |k|i < A: - 1}, 

and the two components are closed and disjoint subspaces of %. 

Step 4 (£-invariance, a{C\'Hk) ): The £-invariance of the finite dimensional com¬ 
bination of eigenfunctions span{^k : |k|i < k — 1} is evident. For every A G cr(£) 
also the corresponding kernel kern£ a is /1-invariant. Therefore the expression 
(|4.26p has to be /1-invariant, since it is just a (finite) direct sum of /1-invariant 
spaces. 

Concerning the spectrum of L in Hk we recall that cr(£|kernz;,x) = o'(£) \ {A}. 
Thus, we obtain from (14.261) that a{C\'Hk) = {—c • k : |k|i > fc}. □ 

After having established the subspaces "Hfc we now turn to the semigroup which 
is generated by £. 

Lemma 4.15. The Fokker-Planck operator C generates a Co-semigroup of bounded 
operators in TL, which is denoted by (e‘'^)t>o. 

Proof. From (14.141) in the proof of Lemma [4.71 we find that for f = i(l-|-,5^-|- 
TrC) the operator (£ — C)lcg°(R") and thus £ — C is dissipative. So we may apply 
the Lumer-Phillips Theorem (cf. Theorem 1.4.3 in [33]) which proves that £ — 
C generates a Co-semigroup of contractions, thus £ generates a Co-semigroup of 
bounded operators in "H. □ 

According to equation (1.2) in [30] the semigroup operators e*^ for t > 0 are 
given by 

(4.27) (e*^/)(x) = [ exp(-VQ-V)/(e‘^x-y)dy, 

(47r)C2 det Qt^ 4 ) 

where Q* = (2C)“^(e^**^ — I). We can equivalently use the following representation 
in Fourier space, which is useful for the subsequent analysis. 

Lemma 4.16. For / G "H and t > 0 there holds 

(4.28) = exp ( - e[{2C)-\l - e-^*^)]^) • /(e-‘^0- 

Proof. If t = 0 the identity (14.281) is obviously fulfilled, so we assume t > 0 in the 
following. For f G TL, (|4.27l) is well defined, and we can write it as 

(e‘^/)(x) = (47r)-"/2(detQt)-'/"e*^^(<(.*/)(e‘^x), 
where ^(x) = exp(—/x^Qj^^x). Using the fact that Qt is diagonal we immediately 
obtain that = (det 47rQt)^/^ exp(—With this we can write the Fourier 
transform of (14.271) as 

= (47r)“”/^detQ7^^^e‘'^''‘^ [ ((/> */)(e“^x) exp(-ix • |) dx 

= (47r)-”/2 det * /](e-‘^4) 

= exp ( - e[i2C)-\l - e-2*^)]0/(e-‘^O- 


So (14.271) and (14.281) are equivalent for all / G "H. 


□ 
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In the next step we investigate the long-time behavior of {e*^)t>o on the sub¬ 
spaces Hk- In the subspaces Hk, the analogue of this analysis was presented in 
Theorem I4.11 vil. Its proof was elementary since the eigenfunctions {/ik : k G Ng} 
form an orthogonal basis of H. But in TL the orthogonality of the eigenfunctions is 
lost, which hence requires more technical estimates of the semigroup. For the rest 
of this chapter they will be mostly based on the representation (14.2811 of {e^^)t>o- 

Proposition 4.17. For every A: S No there exists a constant Ck > 0 such that 
there holds 


(4.29) ||e‘^|„JU(„) < Cfee-‘'=^\ Vt > 0, 

where ci is the smallest entry ofc. 

Proof. We fix A: G Ng and take any / G T-Lk- Our aim is to estimate |||e‘^/|||tj. 

Step 1 (pointwise estimates of /): / is analytic on 11/3/2, and since / G 'Hk we get 
due to 114 . 201 ) that f{^) = 0(|||2) as |||2 —0. More precisely, its Taylor expansion 
with remainder in Lagrange form reads for all ^ G 11/3/2: 


/(I) = ^ for some k. G [0,1]. 

|k| = fc 

Lemma [A.41 provides a uniform bound of |V|/| on I1/3//2, for 0 < /?' < /3. Hence 

(4.30) |/(z)|<C|z|§ ll/lt, VZGH/3V2. 

For estimating the semigroup (14.281) in the norm |||. 1^^ we shall need the following 
estimate for each £ G {1,..., n}: For t > 1 we have 

z:=e-‘^(|±i|e,) GH/3V2, V^gM", 

with P' = < p. Hence, (14.301) yields for all | G K": 


(4.31) 


/(e-*^(|±i^e, 


,-tc 


< C 

< C'e“'='=i* 


(«-§-) i: 
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Step 2 (semigroup estimate): For estimating (I4.28P we compute with (14.311) for 


any £ € {1,... ,n} and for t > 1: 

2 

L2(R") “ 


.F[e*^/](|±i|e,) 


exp 


(|±i|e,) [(2C)-i(I-e-2‘^)] 


.P 2 


< 


exp ( - |^[C-i(I - e-2‘C)]|) |/(e-‘c(| ± i|e, 




(4.32) 


< C 

< c 
= c 


R"- 

2\2fe 


-|€l27C 


/(e-“(«±i|e, 


( I )' 

(I)' 


—2kc\ t 


3-|€l27C 


d| 




2k 


d| 


/ 0 \ 2k 


where yc := (1 — 

Summing (14.321) over all £ G {1,... ,n} we conclude: There exists some C > 0 
such that for alH > 1 there holds 


(4.33) 




V/ G . 


But since (e*^)t>o are bounded operators on H, uniformly for 0 < t < 1 (cf. 
Lemma ld.lbp the above estimate (14.331) holds true for all t > 0 with an appropriately 
large constant C > 0. □ 


With this proposition we conclude the proof of Theorem 14.6 


4.4. The perturbed Fokker-Planck operator. Having defined the extension 
of the Fokker-Planck operator C in H we now turn to the investigation of the 
properties of the perturbed operator £ -|- 0. Note that our x-coordinates are such 
that D = I, and C is diagonal, see the discussion in the beginning of Section IT751 
We make the following assumptions on 0: 

(C) Conditions on 0: We assume that 0/ := -d*f for all / G "H, for some function 
-d : M" —?> K. Thereby the convolution kernel d has the following properties: 

(i) The Fourier transform'd can be extended to an analytic function in 0^/2 
(also denoted by ^), and S G L°°(H^/ 2 )- 

(ii) There holds 'd(O) = 0, i.e. -d is massless. 

(hi) The function 

Jo S 

is analytic in and its real part lies in L°°(H^/ 2 )- 

Lemma 4.18. Under the assumptions (C) the operator 0 has the following prop¬ 
erties in H: 

(i) 0 G 

(ii) For every fc G Nq there holds 0: FLk T~Lk+i- 
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Proof. We start by proving |i|. Due to (C)(i) we have for every f G H that 
•^[0/] = is analytic in and since / satisfies (14.81) we find 

sup ||i5/(--|-ib)||i 2 (Kn) < oo. 

|b|i</3/2 

bGR” 


So, according to Proposition 14.41 0 maps H into H. It remains to show it is 
bounded. To this end we use the norm ||| • |||^, see (14.101) . We start with the 
following computation, where £ G {1,..., n}: 






d^ = lim 

f'y^/3/2 Jrti 


W)(|±i6e,)|^dC 


= IH>I'" 


L°°(ft/3/2) 




dt 


Thereby we have used ([u]) in Proposition l4.4l Note that df is the Fourier-transform 
of an element of H, and thus we may evaluate it at the boundary of 17^/2 in the sense 
of L^-functions. We can repeat this estimate for every £ G {1,..., n} and conclude 
from (I4.10|) that 0 is bounded in H with a norm proportional to ||i?||L~(n^/ 2 )- 
Next we show (jn]). According to (14.201) / lies in Hk iff / has a zero of order 
greater or equal to k at the origin. Now due to (C)(ii) i?/ has a zero of order 
greater or equal to fc + 1 at the origin. Since 0 maps H into H (due to Result (0) 
this shows that 0: %k —>■ 'Hk+i- D 


Corollary 4.19. If 0 satisfies (C) then for every k G Nq the space Hk is invariant 
under £ + 0. 


Proof. This is a direct consequence of Proposition ing and Lemma [4. 181 (p l | above. 

□ 


Throughout the rest of this section we always assume that 0 is such that the 
conditions (C) are satisfied in H for some /3 > 0. Now we fix this /3 and consider 
H with the corresponding weight function w(x) = X]r=i cosh ftxi. In the following 
we discuss properties of £ + 0 in H, which then lead to the final theorem. 

Lemma 4.20. The spectrum a{C + 0) consists entirely of isolated eigenvalues. 

Proof. According to Theorem l4.61 £ generates a Co-semigroup of bounded operators 
in H and has a compact resolvent. Due to Lemma 14.181 (III). 0 is a bounded opera¬ 
tor. Thus we can apply Proposition III. 1.12 in m, which proves that is 

compact for every f G p{C -I- 0). 

It now remains to apply Theorem III.6.29 in [26], which proves that cr(£ -f 0) 
consists entirely of isolated eigenvalues. □ 

In order to characterize the spectrum of £ -I- 0 and the corresponding semigroup 
we introduce the operator dt: H ^ H: f f * 4>. Thereby if is defined by 

:= exp J ^-d(^^s^) ds) . 

As we shall see below, provides a similarity transformation between the resolvents 
of £ and £ -I- 0. 
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Lemma 4.21. 4' satisfies the following properties in 'H: 

(i) For every fc S Nq the operator 4' zs a bijection from Hk to Hk- 
(ii) Both 4' and its inverse 4'“^ are bounded. Thereby 4'“^/ = iF~^[f/ip] for 
all f 

Proof. For the moment we define the operator := iF~^[f/ij}] for all f GFL, and 
show in the following that it is the inverse of To begin with we note that, due 
to the condition (C)(iii), both ip and l/ip are analytic and uniformly bounded in 
^/ 3 / 2 - Thus it follows analogously to the proof of Lemma 14.181 (P that both 4* and 
T are bounded operators in H. 

Since ip and l/ip both do not have any zeros in FLp/ 2 ^ it follows from the char¬ 
acterization (|4.20l) of the space FLk that 4* and th map FLk into itself for every 
k€No. 

Finally we observe that for every f G FL there holds 4'4'/ = 'h'l'/ = /, which 
finally proves that ’F = □ 

Proposition 4.22. There holds 
(i) a{C -I- 0) = cr(£). 

(ii) For every k S Ng the function /k := ^'/ik is an eigenfunction of C + Q to 
the eigenvalue — c • k. Furthermore, for every f G a{L + Q) 

ker(C -{C + 0)) = span{/k : -c • k = C}- 
(Hi) The eigenfunctions /k satisfy /k = V’^/o for all k G Ng ■ 

Proof. Due to Lemma 14.201 we know that the spectrum of £ -I- 0 consists entirely 
of eigenvalues. So, in order to determine the spectrum we look for G C and non¬ 
trivial solutions / G "H of (C — £ — 0)/ = 0. After applying the Fourier transform 
this equation reads 

{c + mi)f+fcyj = df. 

We now make the (non-restrictive) ansatz / = pip. Note that due to (C)(iii) and 
Ip ^ 0 in 12^/2, the requirement f G FL implies that p is analytic in 12^/2- A short 
calculation shows that ip'd = ^CV^ip. Using this, we obtain the following equation 
for p\ 

(C + I4l2)p + I^cvp = 0. 

We find that this is exactly equation (14.151) . In the proof of Lemma 14.81 we have 
shown that 0 ^ p G H is a solution iff C G {~c • k : k G Ng }. And for a fixed C G C, 
p G span{/ik : —c • k = C}. □ 

Note that fo(0) = ip(0)pi.o(0) = 1, hence /o has mass one. 

Proposition 4.23. £-1-0 generates a Co-semigroup of hounded operators, (e**^'^'''®^)t>g. 
For every fc G Ng the space FLk is invariant under the semigroup, and there exists 
some Ck > 0 such that 

Vt>0. 

Proof. According to Proposition 13221 the eigenfunctions of £ and £ -I- 0 are related 
by /k = 'kpk, for every k G Ng. So we find for every ( ^ cr(£) and k G Ng that the 
resolvents satisfy 

FtciO^k = ^ ^ ^ j^ .fk = 
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Since span{/4k : k G Nq} C "H is dense and all operators in the above formula are 
bounded, we conclude the following operator equality in %■. 

(4.34) = i?£+e(C). 


Take any /c G Nq. According to Corollary 14.191 and Lemma [4.211 the identity 
(14.341) holds also in Jik, and Rc+e{C,) is a bounded operator in 'Hk- Now we apply 
the Hille-Yosida Theorem to the decay estimate for (e*'^)t>o stated in Theorem 14.61 
(jyj). It shows that for all m G No and Re^ > —kci there holds 

where > 0 is the same constant as in (14.121) . Applying this resolvent estimate 
to (14.341) yields for all m G Nq and ReC > — fcci: 


\\Rc+e{Cr\nJ 


nk\\m(Hk) ^ 




(ReC + A:ci)'" 

Applying the Hille-Yosida Theorem again implies that £ -I- 0 generates a Co¬ 
semigroup of bounded operators, which satisfies the following estimate: 


where 0 < C* < Cfe||4'||^(«,)||4' ^Wss^Uk)- 


□ 


We conclude this section by summarizing the main results. 

Theorem 4.24. Under the conditions (C) on 0, the perturbed Fokker-Planck op¬ 
erator £ -G 0 has the following properties in 'H: 

(i) a{C -G 0) = cr(£) = {—c • k : k G Nq}, i.e. £ -G 0 is an isospectral 
deformation of C. 

(ii) The functions /k := di/ik are eigenfunctions o/£ -G 0 for all k G Nq. For 
every A G cr(£ -G 0) the corresponding eigenspace is given by 

ker(A — (£ -G 0)) = span{/k : —c • k = A}. 

(in) For every k G No, the operator C+Q generates a Co-semigroup (e‘*^'^+®^)t>o 
on TLk, and there exists some constant Ck > 0 such that 

w>0^ 

In particular, this theorem implies exponential convergence of the solutions of 
the perturbed Fokker-Planck equation towards the stationary solution: 

Corollary 4.25. Let ip he given, and let f{t) := he the corresponding 

solution of dSD. Set m := dx G C. Then there exists a constant C > 0 

such that 

11/(0 - mfoWkj < C||(^- Vt > 0, 

i.e. f(t) converges exponentially to mfo. 


Proof. Since /o is the unique normalized zero eigenfunction of £ -G 0 we obtain: 

fit) - mfo = e*(^+®0£ - mfo). 

Since (f — mfo has zero mean, it follows from Lemma 14.121 that it lies in 'Hi. But 
(e‘('^+0))t>o decays exponentially on Hi with rate —ci, see Theorem 14.241 (Iml) . So 
we get for alH > 0: 

11/(0 - mfoWu, = ||e*(^+0)(v?- to/o)|U < CiIKg^ - m/o)||^ e“‘'=L 
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□ 

Remark 4.26. Note that £ + 0 is neither self-adjoint in H nor in TL. But the fact 
that a{£ + <d) C M and that £-)-0 is only a “deformation” of £, see (I4.34L suggests 
that £ -|- 0 is self-adjoint in an appropriate space. To verify this we introduce the 
inner product 

jR" M 

and the corresponding norm || • Hy,. The associated space Sj is the set of all functions 
such that II • ||i 5 is finite. This is indeed a Hilbert space, and 4* is an isometry between 
H and Using (I4.34p we see the self-adjointness of L -I- 0 in ij: 

{{L + e)f,g)s^ = (ToLo4'-i/,5)f> 

= (/) {L + Q)g)sj, 

where we have used the self-adjointness of L in H. In the eigenfunctions /k of 
£ -I- 0 are orthogonal again (like the functions /Xk in H). Altogether, we conclude 
that L in H and £ -|- 0 in .5 are isometrically equivalent via the map 4/. Hence, 
£ -|- 0 inherits most properties of £. However, we point out that discovering the 
map 'k, without the preceding analysis, is a non-trivial issue. 

Furthermore, the Hilbert space is difficult to be characterized explicitly. In 
particular, it is usually not possible to describe Jo as a weighted £^-space. A 
simple calculation shows that = L^{v) for some weight function v only if for all 
/ G C^CR'^) there holds 

But in general this function v will not be independent of /. 


Appendix A. Results in functional analysis and deferred 

proofs 

On fl = R” it is possible to find compact embeddings of weighted Sobolev spaces 
into weighted £^-spaces if certain conditions on the weight functions are satisfied. 
Here we need the following corollary from Theorem 2.4 in [35]: 


Lemma A.l. Let v,w be two weight functions on LI = R”. Assume further that 

,, , iy(x) 

(A.l) hm ess sup = 0. 


Then there holds the compact embedding H^{v,w) Lp‘{w). 

Lemma A.2. Let v be a weight function on MT. Then (7(j”(R") is dense in Lf{v). 

The proof of the above lemma is straightforward, see [36] for more details. 
Lemma A.3. There exists a constant C > 0 such that for every f € TL we have 

|VV(0)I < CII/IU, kGN^. 
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Proof. We have 

= ||^[xV(x)]|U~(Rn) < ||xV(x)||il(Rn) 

= / l/(x)|w(x)^/2 • |x’"w(x)"^/2|dx < ||/||<^( f x2*"a;(x)"Mx) ' . 

Since a;(x)“^ decays exponentially as |x|2 —>■ c» the last integral on the right hand 
side is finite. □ 


Lemma A.4. For every 0 < f3' < /3 and k G Ng, there exists a positive constant 
C such that 

sup |vV(z)|<c||/iu, yfen. 


Proof. Due to Proposition l4.41 all functions f G PL satisfy /(4+ib) = J^[/(x)e’^ ’‘](|) 
for |b|i < /3/2. Hence, 

(VV)(4+ ib) = J-[(-ix)‘^/(x)](^+ ib) = J-[(-ix)V(x)e^-](0 


follows for |b|i < ff /2 and k G Ng . Then, 

sup |V*"/(z)| < sup ||J'[(-ix)'"/(x)e'"'’"](|)||i=o(Rn) 

zeap,/2 \h\i<p'/2 

< sup ||x‘"/(x)e'^'’"|lLi(Kj) 

|b|i</3'/2 


The norm " 




can be estimated as 


^2kg2bx 


V a;(x) 


dx < 


f x 2 ke/ 3 'I^U 
V w(x) 


dx < 2 




(fx) 


;(x) 


dx =: < oo, 


where C is finite due to 0 < /3' < /3. Thus, the estimate sup^gQ^,^^ |V*^/(z)| < 


C WfWui for all / G 'H follows. 


□ 


Lemma A.5. Let X ^ X be Hilbert spaces, and i/jq, ..., i/'is-i G .^(A,C) be 
linearly independent functionals. Then ipj := C) for allO < j < k—1, 

and 

k-l k-1 

Pi ker -ipj = 01^^ ker -ipj. 
t=o i=o 

This result coincides with Lemma C.2 in m- The proof can be found in therein. 

Lemma A.6. Consider two Hilbert spaces X ^ X and a projection Px G 
such that Px := Pa|x G L^{X). Then ranP^ = cl^ ranPx and kerP^ = cl^ kerPx. 

This result coincides with Lemma C.l in m- 
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